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1. SUMMARY 


This volume of the report describes the data processing procedures and 
the computer programs developed to predict structural responses using the 
Impulse Transfer Function (ITF) method. There are three major steps in the 
process : 

• Analog- to- digital (A-D) conversion of the test data to produce 
Phase I digital tapes 

• Processing of the Phase I digital tapes to extract ITF's and 
storing them In a permanent data hank 

• Predicting structural responses to a set of applied loads 

The analog to digital conversion is performed by a standard package which 
will be described later in terms of the contents of the resulting Phase I 
digital tape. 

Two separate computer programs have been developed to perform the 
digital processing: 

Program I 

• ITF Program - extracts ITF’s and stores them in the data bank 
Program II 

• Response Program - predicts structural responses to a set of input 
forcing functions 

All coding was initially done in Fortran IV for the IBM 360/T5; this 
program was used to produce the results presented in Volume I. The data bank 
used was an IBM 231^ Disk Pack, with random access capability. Both programs 
have been modified to Fortran V for the Uni vac 1108, using standard 7-track 
tape as the data bank. 

The programs were originally sized for the present problem of 36 applied 
impulses and 70 response channels. The tape storage version, however, allows 
a virtually unlimited number of responses and the number of applied impulses 
can be increased with only minor program revisions. 
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2. BACKGROUND AND INTRODUCTION 


Analog Tape 

Test data was recorded on 108 KHz vi deband analog tape at 60 lps. 

This tape has l4 channels, of which one is used exclusively for an IRIG-B 
time code, one for the 108 KHz reference, and another for the impulse. The 
ITF's are multiplexed, 10 to a channel, allowing a maximum of 110 ITF's per 
tape. For the LTA-11 testing, TO responses were recorded on channels 2-8, 
the impulse on channels 1 and 12, IRIG-B on channel 14, and the 108 KHz 
reference on channel 13. Channels 9 and 10 were not used. Time sequencing 
of events as recorded on the analog tape is schematically shown on Figure 2 
and described in the accompanying table. 

Time (Seconds) 

0-9 
9 - l4 
Ik - 19 
19 - 24 

24.00 - 24.01 

24.01 - 24.200 
24.200 - 30. + 

Oscillograph playbacks of each ITF was made to verify the data and for 
later comparison with plots made of the digitized data. For the present con- 
tract, 36 impulses were applied to the vehicle; one analog tape was produced 
for each applied impulse. Sixty-eight active measurements per impulse were 
recorded for a total of 2448 ITF’s. 

An alog-t o-Digital Conversion 

Conversion of the analog tapes was performed by the LM Data Reduction 
Station (LDRS) using a general A-D conversion program. Seven Phase I digital 
tapes are produced from one analog tape, with each containing the impulse, 
and 10 ITF's. Conversion starts at the beginning of the D.C. calibrations 
and continues for at least 6 seconds. Most of the LTA-11 conversions were 


Signal 

A.C. calibration signal 
Hard zero, gauges shorted 
D.C. calibration - 100? of full scale 
Hard zero, gauges shorted 

10 ms event pulse - approximately 100? full scale 
Recorded zero, gauges in circuit 
Impulse occurs 200 ms after event pulse and data 
is recorded for at least 6 seconds 



A.C. CALIBRATION 


D.C. CALIBRATION 


1 sec 



FIGURE 2 

ANALOG TAPE Schematic - Response Track 
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perforated for only 3 seconds due to the storage limitation in the earlier 
version of the computer system. Subsequent modifications have increased 
storage to 6 seconds. During conversion, the impulse track is sampled at 
20,000 samples /sec, and the ITF's at 5000 samples /sec. Formatting of the 
Phase I digital tape, which must be compatible with the IT? program (Program 
I), is discussed next. 

Phase I Digital Tapes 

Phase I digital tapes produced by the A-D conversion are standard 7 
track input tapes for the ITF program. Each tape contains the time history 
of the impulse track identified by the load point, and ten of the 68 asso- 
ciated ITF's. The histories were directly digitized from the analog records 
starting with the D.C. calibrations and continued for 3 seconds. In order to 
satisfy the requirements of the ITF program the following specifications must 
be adhered to for the Phase I digital tape: 

• It must be a 7 track binary tape containing two files of data, 
separated by a standard End of File mark. 

• File 1 contains addressing and record size information which is 
not used by the ITF program. This is generated by the generalized 
A-D conversion program and is passed over by specifying on the 1/0 
control cards that the second file is to be used as the input data 
set. 

• File 2 contains the time history data written in records, each 
4868 characters long. This may be interpreted as: 

29,208 bitB 
9*736 octal digits 
4,868 BCD characters 

2,434 12 bit data values (4 octal digits) 

• The data is most meaningfully interpreted as 12 bit words ranging 
from 0 to 7777 octal (4095 decimal). 

• Each record represents 32 milliseconds of real time. 
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Data 

Status and time values from IRIG— B 
640 values of Impulse track (32 ms x 20 samples/ms) 
160 values of response 1 (32 ms x 5 samples/ms) 

160 values of response 2 
160 values each of responses 3-10 
Dummy values to fill out record 

• Each record contains data for 10 response channels; an eleventh 
response channel could be accounted for but its inclusion would 
cause some added complexity in the numbering scheme* 

• Data should start (Record l) within ^3*2 seconds (100 records) of 
the start of the DC calibrations (i.e., no more than 100 records of 
zeroes at start, and at least 50 records of DC calibrations) 

• There should be at least 75 records (2.4 seconds) of zeroes between 
DC calibrations and event pulse. Normally, with direct conversion, 
there is a period of 5 seconds (^ 156 records) between DC calibrations 
and event pulse. In the present contract work, some conversions were 
performed with time Jumps in this zero area. This is not considered 
to be a useful economy, but can be accepted by the I TP program. 

e Data should be recorded continuously from the event pulse to the end 
of data. For the LTA-11 only 3 seconds of response was recorded, but 
present storage capability for the 1108 program is 6 seconds. 

• During A—D conversion, the DC calibrations are used as an indication 
of the full range of values required. Scaling should be such that 
the DC calibrations fall at approximately 90* of the full digitized 
scale. 

• No sign bit is used, hence the digital value of 0 is equivalent to 

- 100 *. 


Decimal 

Octal 

Reading 

0 

0 

-100* 

2048 

4ooo 

0 

4095 

7777 

+99.95* 


• Scale factors of all measurements must be supplied to the ITF program 
to allow eventual conversion to engineering units. 
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• DC calibrations, the event pulse and the impulse are to be recorded 
in the positive direction during recording of the analog tape and 
for A-D conversion . Sign changes to account for direction of the 
impulse relative to the coordinate system are data input to Program IX. 


Program I - ITF Processing 

This program extracts the ITF^ from the Phase I tapes, evaluates the 
zero adjustment and scale factor, locates time zero, and stores the ITF f s in 
the data bank. Plotted output of the ITF's is available for visual inspection 
and comparison to the oscillograph records; conventional computer printout is 


also available. 


How the scale factors and zero adjustment are obtained is discussed next 
with the aid of a specific example. (See Figures B-2, -3, -4) 

Consider response 1, vhich is gauge GA-12. It is a force gauge located 
on the aft interstage truss. It is calibrated prior to testing, and it is 
deterained that its output, when passed through a Preston Filter, with gain 
set at 1000 and an RC of 100K, is 2.5 volts D.C. at 5360 #. The nominal 
range of the analog recorder is +2.5 vdc. From a sample run it is observed 
that the output of GA-12 is low; the RC of the Preston Filter is then reset 
to 4J0K, at which value the output of GA-12 is 2.5 vdc at 1143 Since 
this output signal is satisfactory, it is duly recorded that 1143 i? a 1003» 
full scale * 2.5 vdc. This data defines the engineering units of the measure- 
ment. Following an application of an AC calibration signal to the recorder, 
GA-12 is in the circuit and recording at a nominal zero level for 5 seconds, 
then a direct input of 2.5 vdc is applied for 5 seconds, followed by 5 seconds 
of a "hard 11 zero level during vhich the gauge is shorted out of the circuit 
(see Figure 2). An event pulse of 2.5 vdc is then applied for 10 msec and the 
impulse and recorded response occurs 200 msec later. During the period between 
the event pulse and the impulse, gauge GA-12 is in the circuit and its output 
is being recorded. Examination of the Phase I digital tape in this section of 
the test shows the following (adjusting counts to the range -2048 to +2047)s 
Test Period Record Reading (average) 

End of DC cals 250 - 300 1617 counts 

Start of hard zero 310 - 350 8 counts (gauge shorted out) 


Event pulse 463 513 counts 

Event to impulse 464 - 468 10 counts (gauge in circuit) 




8 


From this we can surmise that: 

1143 0 = 2.5 vdc = 1609 counts (l6lT-8). 

It is important to note how the hard zero and applied 2.5 vdc pulse are used 
to obtain the scaling. 

GA-12 is recorded at 10 counts, while in the circuit, prior to the 
application of the impulse. This reading represents 1.4 0 , and is computed 
from the tabulated values . 

(10 - 8) x HU3/1609 * 1.4 

This indicates that the strain gauge, within the supported structure, 
prior to impulse application is not in a wholly unloaded state. This level 
of 10 counts is taken as the nominal zero level of the measurement (zero 
adjustment). Further complications concerning oscillations of the zero level 
due to swaying of the test article will be discussed later. To complete 
these observations, following the impulse, GA-12 records oscillations which 
range between +183 counts and -209 counts and eventually settle down to a 
small range around 10 counts. These peak values represent: 

(183 - 10) x 1143/1609 - 123 0 
(-209 - 10) x 1143/1609 * -155 0 

This procedure is also followed on the impulse track, to obtain the load 
time history of the impulse. In a later calculation, this data is used to 
compute the impulse magnitude (area under the load time history) which is 
eventually used to scale the impulse response to a unit impulse. 

With the above discussion in mind, the processing logic of the ITF 
program may now be considered. Figure C-2 shows a flow diagram of the sequence 
of operations in which one Phase I tape, containing an impulse track and ten 
responses, is processed to obtain the ten ITF f s. Records are processed sequen- 
tially. Note that each record contains: 

• 32 milliseconds of real time 

• 640 impulse track values (20 samples/ms x 32 ms) 

• 160 values per ITF (5 samples/ms x 32 ms) 

The programmed procedure is as follows: 

1. Nine records are skipped, the tenth is read and the average values 
of the impulse and response tracks are calculated and saved. In 
reading records, counts are automatically adjusted to the range 
-2048 to +2047 (= 100£ full scale). 
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2. Step 1 is repeated until a record is found where the impulse track 
average value exceeds 1000 counts. This is the region of DC calibration. 

3. Every tenth record is averaged and the averages saved, until a record 
is found whose impulse track average is less than 500 counts. This is 
the hard zero region following DC calibration. 

4. Fifty more records are passed, with every tenth being averaged and 
saved. The last of these is used to establish the zero level for 
locating the event pulse. 

5. One-hundred fifty records are averaged and saved. Each is examined 
to locate one where the average of the impulse exceeds the zero 
established in step ^ by 120 counts. This locates the record 
containing all or most of the event pulse. This record is checked 
to determine that the sum of three consecutive values exceeds 2000. 

6. Seven records are read, averaged and the averages saved. The last 
two records are saved in their entirety since they contain the 
impulse and start of the ITF. 

7. The two records containing the impulse are examined further. Using 
the average of the three previous records as a relative zero, three 
consecutive values over 200 counts are sought to indicate the start 
of the impulse. 

8. Time zero for the ITFs is taken as the nearest millisecond to the 
first 200 count value of the impulse. 

9. The area of the impulse is computed using the preceding 2 milli- 
seconds as the relative zero level. 

10. With the time zero established, values from the ITF tracks are 
stored from the two records previously read. 

11. Additional ITF records are read and the values stored until 6 
seconds of response is obtained or the end of the data reached. 

12. Values are stored as 1 ms averages, with every 5 points being 
averaged. Values are stored without adjustment to the zero level 
or the scale factor. 

13. Scale factors for the impulse and ITF’s are calculated using the 
average record values saved from the 20th, 30th and Uoth records 
after the start of DC calibrations and the 10th, 20th and 30th 
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records after the end of DC calibrations . Engineering unit values 
of the impulse and ITF tracks corresponding to 2.5 vdc (100# scale) 
are supplied as input to the program. Coupling the engineering 
unit value, with the area under the pulse, the final scale factor 
is computed in engineering units per count for a unit pulse. 

14. This scale factor is now stored in a status file along with other 
information such as ITF duration, impulse area, DC calibration step 
size, engineering units full scale value, and two zero values, one 
calculated from the average of the records between the event pulse 
and the impulse, and the other from the record averages of the ITF’s. 

15. Printed output includes pertinent record averages throughout the 
process as well as several records in their entirety. 

ITF Storage 

The ITF processing subroutine assumes that the ITF’s are stored in tempo- 
rary, direct access data files on Fastran Drum, referenced by Define File state- 
ments in Fortran V under the EXEC -8 operating system. The use of disk pack or 
tape for permanent storage Is discussed later. 

A full length ITF Is 6000 values in 1 ms increments. ITF values are read 
from the Phase I tapes as integer values in the range of 0 to 1*095* They are 
adjusted to the range -2048 to +2047 for ease of analysis and interpretation. 
Five point averages are computed; these are also in the range of -2048 to +2047* 
To store the 6000 values economically, they are reconverted to the range 0 to 
4095> which is represented by exactly 12 bits, and packed three to a 36 bit 
word before storing as a 2000 word string. The file for ITF storage consists 
of 3600 - 200 word records. One ITF is ten records. The file holds 360 ITF’s, 
representing ten ITF’s for each of the 36 applied pulses. Ordering of the ITF’s 
is by measurement number, l.e.» ITF 1, Pulses 1-36; ITF 2, Pulses 1-36, etc. 

The values stored are not adjusted for the nominal zero and are not scaled to 
engineering units. The zero level and scale factors are stored in a second 
file of 10 - 216 word records. This is, in effect 360 - 6 word arrays, ordered 
by ITF number and loads the same as the ITF’s. These six words contain float- 
ing point values of ITF duration in seconds, two zero values, the scale factor 
from counts to engineering units, impulse area, and the impulse scale factor. 

For computing economy, the ITF’s are scaled in Program II when predicting 
responses * 
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Program II - Response Program 

The response to a set of forcing functions is obtained through the appli- 
cation of the Duhamel integral, as described in Volume I. The numerical solution 
proceeds as follows: 

• An applied load time history is read as input. 

• A time interval is specified as input; its selection depends on the 
expected frequency content of the output. 

• Total time of the response is specified as input. 

• The applied force time history is interpolated to each time interval 
and saved as an array. 

• The ITF resulting from the application of a specific impulse is read 
in from the data bank along with zero adjust and scaling information. 

• The ITP is passed to the filtering routine as specified by the user. 
(Filtering is discussed in a later paragraph.) 

• The filtered ITF is reduced to a set of averaged points corresponding 
to the time increment over the total time. 

• Zero adjustment is made. 

• The response is calculated from: 

k 

R(t k ) = cat E Wi - i 

i=l 

where for a typical gauge reading force, 

R is in lbs 
At is in sec 

F is in lbs (forcing function) 
h is in counts (the ITF) 

C is in lbs /lb -sec- count , (scale factor) 

• This process is repeated for all other load points, with options to 
plot and print individual and/or summed responses. 

Filtering 

ITF filtering is not an integral part of the ITF method, but vas intro- 
duced to eliminate spurious low frequencies in the measured ITF's arising from 
the LTA-11 supports. These effects are removed from the ITF's during response 



12 


processing by a filtering routine utilizing a Fast Fourier Transform technique, 
developed at Grumman. 

During performance of the present contract work, these low frequency 
oscillations were observed in the ITF’s and several methods were considered 
to eliminate them. The most useful was the Fourier Filter technique which 
did not require visual inspection of the ITF plots and was similarly applicable 
to the different types of responses (forces, strains, accelerations). 

Filtering is done in the response program following retrieval of the 
ITF’s from the data bank, and prior to their use in calculating a response. 

The ITF is passed to the filtering routine where the following opera- 
tions are performed. 

• Stored data that defines an ITF time history (maximum = 6000) are 
time averaged to reduce the number of values to a more workable 
number (typically 300). This averaging eliminates the higher 
frequencies, but retains the character of the low frequencies 

of interest here. 

• The reduced number of time history points are transformed to the 
frequency domain. 

• All Fourier coefficients above a specified cutoff frequency are 
zeroed, and the points transformed to the time domain. They now 
represent the low frequency content of the ITF as a time history. 

• This low frequency time history is then interpolated back to the 
original sample rate and subtracted from the original ITF to obtain 
a new function which is considered filtered. 

• When filtering, it is not necessary to zero adjust the ITF, since 
the process removes any biasing in the data. 

Zero Adjustment 

The need for accurate evaluation of the zero level of the ITF’s can be 
recognized if it is recalled that responses are obtained by integrating the 
ITF’s. An error in the zero level will result in a steadily growing error 
over the duration of the response. Three methods of zero level evaluation 
are available, under input control. 
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• When processing the ITF from Phase I tapes, the calibrated zero 
level, which exists between the event pulse and the time of impulse 
application, is determined and stored in the status file* This 
zero level is of short time duration and may not be sufficiently 
accurate. 

• A second estimate of the zero level is obtained by averaging the 
ITF over all but the early time period. The time skipped should 
be equal to several periods of the lowest natural frequency of the 
test article. This zero level will assure that the response to a 
constant load does not change its gross level significantly after 
the first few low frequency cycles . This estimate of zero level is 
done in Program I and is stored in the status file.- 

• The third method of zero adjusting is performed automatically whenever 
the filtering routine is used. This routine removes all frequencies 
below a specified value, including the zero frequency. 

Data Bank 

The system was originally designed using an IBM 23lk disk pack, capable 
of permanently filing 2520 - 80,000 bit ITF's. Since this capacity for per- 
manent storage does not at present exist on a direct access data device at 
the MSC Uni vac 1108 facility, the system was redesigned to use tape storage. 
The primary progranming problem, in such a system, is to minimize the number 
of tape operations for efficient use of computer time. This was done by 
making use of a smaller direct access temporary data set on Fastran Drum. 

All programs, and the associated data files have been sized to a 10 response 
system, with provisions in each program to initialize the data files from 
tape, and, at the end of an operation with any set of 10 responses, to dump 
the files back to tape. It is necessary to work with 10 responses simul- 
taneously because 10 responses are grouped on each Phase I tape. When exe- 
cuting the ITF processing program, if severs! Phase I tapes for the same 
set of 10 responses are processed sequentially, drum to tape 1/0 does not 
consume a major portion of operating time. 

When executing the response program it is only necessary to load the 
drum from tape. Final dumping is not required, since the permanent data 
base is not modified. All subroutines operate on the basis of ten numbered 



responses, while the main programs, which perform the loading, relate the 
responses to the overall numbering scheme . If large scale permanent direct 
access storage becomes available, either the tapes can be replaced with the 
permanent data set and the programs operated exactly as with the tapes, or 
else the temporary data sets may be increased to full size, with minor pro- 
gram modifications. 
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3a Program Description - ITF Program, Program 1 

Program I, the ITP Program, processes Phase I digital tapes to extract 
the ITF’s, stores them in the data hank, and performs other data bank main- 
tenance and utilization functions. 

Figure C-l shows a simplified flow diagram of the main program. It 
performs CLEAR, LOAD, DUMP, STATUS and PLOT operations on the data bank, 
and calls subroutine STORED, which contains the bulk of the programming, to 
process Phase I tapes. To understand these functions it is necessary to 
review how the program handles the data bank described in Section 2. 

Program I references two temporary, direct access data files through 
define file statements. 

DEFINE FILE 1 (10, 2l6, U, IUl) 

DEFINE FILE 2 (3600, 200, U, IU2) 

File 1 contains data describing the status of 360 ITF’s. The ten records 
pertain to 10 responses. The 2l6 words per record are an array of (6, 36) 
real variables, six words for each of 36 load points. 

word 1 - conversion factor, units per count, which converts the value 

of the stored ITF from digital counts to the engineering units 
associated with the particular response point 
word 2 - the zero level of the ITF, derived by averaging the analog 

response for the period between the event pulse and in impulse 
word 3 - the duration of the stored ITF in milliseconds; maximum ** 6000 
word h - the zero level of the ITF derived by averaging the ITF from 
320 milliseconds to the end 

word 5 - scale factor for the impulse track in pounds per count 
word 6 - area of the impulse in pound- seconds 
File 2 contains the 360 ITF’s, ten records per 2000 word ITF. They are 
ordered by response point (i.e.. Response 1, loads 1-36; Response 2, loads 
1-36, .... Response 10, loads 1-36). 



For a specific problem with 10 responses or less, permanent data files 
could be used and the processes of dumping to, and loading from tape could 
be dispensed vith. The programs can, however, deed with an unlimited number 
of response points, by the use of binary tapes to store the contents of data 
files 1 and 2. Through use of the LOAD and DUMP operations, the contents of 
these files may be dumped to a binary tape or loaded from a binary tape, with 
no real restrictions on the number of tapes being used for a problem. 

The effects of the five data bank operations are as follows . 

CLEAR - The contents of files 1 and 2 are set to zero. This operation 
is usually performed when the first Phase I tapes for a set of 
responses are being processed and no LOAD operation is performed* 
LOAD - Contents of a binary data bank tape are loaded into files 1 and 2. 
Presumably, this tape was created by a dump following an earlier 
run. 

The second content of the binary data bank tape is: 

record 1 - 100 words - an integer variable representing the 

tape sequence number - repeated 100 times. 

record 2 - 2160 words - the entire contents of file 1. 

records 3 - 362 - 2000 words each - the entire contents of file 

2, but in 2000 word records rather than 200. 

DUMP - Contents of files 1 and 2 are dumped to binary tape to create a 
new or updated data bank tape 
STATUS - prints the contents of file 1 
PLOT - plots or prints specified ITF’s, stored in file 2. 

The sequence number on the binary tapes is used to insure that the proper 
binary tape is mounted, since a number of them will be used for the typical 
Job. While most of the internal logic of the program deals vith ten responses 
numbered 1 to 10, an overall numbering system should be employed, whereby the 
responses are numbered consecutively from 1 and are stored vith tape number 1 
containing responses 1-10, tape 2 containing 11-20, etc. Input data pertain- 
ing to response numbers always refer to the overall numbering scheme. This 
is converted by the program to a tape sequence number and the internal number 
range 1-10. 
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An additional storage operation vas written into the program in order to 
transfer data between Grumman *s disk pack system and MSC’s tape storage system. 
This operation is: 

LOADGQ - loads an existing BCD tape containing a number of ITF's. 

The END operation causes a normal termination of the program. 

Processing of Phase I Tapes 

The RUN operation is used to process Phase I tapes. This function is 
performed primarily in subroutine STORED which also calls RITREC, READ7» AVG 
and UNPACK. Figure C-2 shows a flow diagram for subroutine STORED. The 
basic sequence of operations in processing the Phase I tape is covered in 
Section 2. As each record of 32 milliseconds is read, it is separated into 
impulse track values and 10 sets of response track values. The average value 
of each of these sets is determined for each record examined. Only these 
averaged values are used until the actual impulse is encountered and the 
responses start . 

While examining records to locate the DC calibrations and event pulse, 
only the impulse track is examined since events occur simultaneously on each 
track. Once the record containing all or most of the event pulse is located, 
the impulse can be readily found. It occurs 200 ms after the event pulse and 
must be in the 6th or 7th record following the event pulse. The first 5 of 
these are read and their averages saved to use in determining zero levels. 

The 6th and 7th are examined to locate the impulse. 

A temporary zero level for the impulse track is established using the 
average value of the three previous records. The impulse 1 b located by 
finding three consecutive values, each exceeding the average by 200 counts. 

A new zero level is established to use in calculating the impulse area. It 
is based on the 40 values preceding the start of the impulse (2 ms). The 
first and last points on the Impulse whose values exceed 50 are located, and 
a Simpson’s Rule integration performed between them. If the number of inter- 
vals is even, the first interval is evaluated by the Trapezoidal Rule. One 
additional point on each end of the impulse is now used to calculate the areas 
of the end pieces. 



Scale factors for the impulse and response tracks are calculated and 
the impulse scale factor is used to get the impulse area in lb— seconds. It 
should be noted that for this problem, there are 7 phase I tapes for each 
impulse, and that the impulse track on each of them is derived from the same 
analog track* Any differences between them are due to the precision of the 
A-D conversion, and comparison of the impulse areas can serve as a check on 
the A-D process. 

At this point, 64 ms of data, is present in core. The starting point 
of the impulse is known, and this is adjusted to the nearest millisecond, 
to define time zero for the responses. The remaining task is to transfer 
data from the response tracks of the Phase I tape to the data file in which 
the ITF*s are stored. To do this, the intermediate array IDISK is established 
dimensioned to (600, 10). It holds 600 ms of ITF data for the 10 responses 
on the phase I tape • Records are now read into core through the use of the 
JREC array, converted by averaging from 5 points /ms to one point/ms and 
transferred to the IDISK array. The process is started by transferring the 
appropriate values from the two records containing the impulse. When the 
IDISK array is filled, it is packed down to (200, 10) since ITF's are stored 
as three values to a word. The 200— word columns are then stored in the 
appropriate records of file 2 and the process repeated. Processing continues 
until 6000 milliseconds of data is stored, or the data runs out. A record of 
the number of milliseconds of valid data is kept for inclusion in the STATUS 
file (File l). When processing of the phase I tape is completed, the status 
items are loaded into File 1, summary output is written, the tape is rewound, 
and control is returned to the main program. 

A summary of the subroutine cadis follows. 

SUBROUTINE STORED (ITAPE, IRESP, IRESPP, ITEST, FULLI, FULL, NG) 

ITAPE - fortran unit number on which phase I tape is mounted 
IRESP, = range of response numbers to be stored, inclusive. They 
IRESPP are internal response numbers in the range 1 to 10 
ITEST » a value (from 1 to 5) that controls the number of records 
written as printed output 
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HJLLI a scale factor in pounds per full scale for the impulse track 
FULL = 10 item array with the scale factors in engineering units per 
full scale for the 10 response tracks 
NG » a counter set to zero after successful completion 

SUBROUTINE RITREC (JB, IB, IC) 

JB » record number 

IB * 2400 word array containing the impulse and response values 
from the record adjusted to the range of -2048 to +2047 
IC = (l60, 11 ) array containing the responses. 

This routine prints the contents of one phase I tape record. 

SUBROUTINE READ! ( JTAPE, IREC1, IREC2, IAVG, KS, NG) 

JTAPE = fortran unit number 

IREC1 « 34 member array - not used 

IREC2 = 2400 member array 

IAVG « 11 by 500 array 

KS = record number 

NG = counter 

This routine is used to read one record from the phase I tape. The 
actual read is performed by the MSC library routine MREAD, used to read non- 
standard tapes. MREAD is called in UNPACK which is called by READ?. The 640 
word set of impulse values and 11 to l60 word sets of response values sure 
loaded into IKEC2. 

Further processing is under control of the input, NG. 

NG * 0 for no further processing 
NG * 1 for impulse track only 
NG = 2 for responses 
NG = 3 for both 

20 48 is subtracted from each word to modify the range of data from (0 to 
4097) to (-2048 to +2047). 

Averages sure performed on the impulse and responses suid stored in the 
IAVG array according to record number. 



SUBROUTINE AVG (IREC, IAVG, NG) 

IREC = 2^00 values of tape record 
IAVG = 11 averages to be calculated 
NG = same as in READ7 

Averaging described in READ7 is performed. 

SUBROUTINE UNPACK (INI, IREC2, ISTAT) 

INI = fortran unit number 

IREC2 * 2^00 values to be returned 

ISTAT * status of read 

As described in Section 2, a phase I tape record contains 2^36 12 bit 
words. These are read by MREAD into an array of 812 - 36 bit words. The 
last 800 of these are unpacked, 12 bits at a time and stored in IREC2. 

SUBROUTINE CRUNCH is only used to decode the special tape used to 
transfer ITF’s to the data bank for demonstration purposes. It may be 
replaced with a dummy routine, or reference to it in the LOADGQ section of 
the main program may be deleted. 

SUBROUTINE PL0T1 ( I FILE, INN) 

I FILE is a work area of 3000 words supplied by the main program 
INN array of 8 control constants 

PL0T1 reads ITF's from data file 2 and either plots them, prints them 
or both. 

Plot arrays X and Y are completely set up, but no calls to plotting 
hardware are supplied. 

Provision is made for calling plot initialization and plot termination 


routines . 
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3b Program Description - Response Program; Program II 

Program II, calculates the structural response to a set of forcing 
functions. It consists of a main program, a filter routine, FILT, a fast 
fourier transform routine, FASTI, and a plot routine, FLOTB. Flow diagrams 
for the main program and the filter routine are shown in Figures C-U and 
C-5. 

The main program uses two data files which are identical to those used 
in Program I. They contain ITF's for 10 response points and 3 6 load points. 
Permanent storage of the ITF’s is on the binary data bank tapes created by 
Program I. The appropriate tape must he read into the data files at the 
start of execution. In order to run program II it is necessary to have pro- 
cessed, for the responses desired, those load points at which forcing func- 
tions are to be applied. Other load point ITF's need not have been processed, 
since the data hank tape contains dummy values where none have yet been 
supplied, so as to keep proper sequencing of data. Input data to the program 
consists of a set of control cards with a function name in columns 1-6 
followed by additional data. 

For the first pass through the program the order of data is fixed, but 
additional passes can be made on one run, with some variation allowed. In 
addition, a set of tables of the forcing functions must be supplied. 

The function names on the control cards, in normal input order are LOAD, 
TIME, FORCE, RESP, PLOT and CYCLE. Forcing function tables follow the force 
card. 

Except for recycling for multiple problems, the flow of the program is 
straight-through, and is best described in relation to the input data. 

The LOAD card specifies an I/O unit number and response numbers and 
results in a data bank tape being read into Files 1 and 2. All features of 
this operation are identical to the LOAD operation in program I. 

The time card specifies the number of points, or time slices required 
by the analysis, and their spacing in milliseconds. The number of points 
should not exceed 1200 and the total time should not exceed 6000 ms. 
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While the ITF f s are stored at 1 ms time intervals, it is not necessary 
to evaluate the structural responses at 1 ms intervals . Selection of the 
calculation interval should he based on system frequencies, hence if the 
highest system frequency is 50 cps (20 ms period), spacing of 4 or 5 ns will 
adequately represent the results. The forcing functions will be evaluated 
at the specified times and the ITF f s averaged to that time interval before 
performing the Duhamel integration. 

The force card is the principle heading card for the forcing functions. 

It is followed by a set of cards describing the forcing function time histories 
at each load point required by the problem. To minimize the effort in pre- 
paring data for the forcing functions, all forcing functions are described 
for the same time points. The time scale is described on the FORCE card as 
a sequence of points at a fixed time interval, followed by another set of 
points at a different time interval, up to 6 sets. Hence a 500 ms time scale 
might be specified as 20 points at 2 ms, 6 points at 10 ms and 4 points at 
100 ms. 

An array is established giving the time value for each point on the 
tabulated forcing function. The force card is followed by a number of sets 
of forces for each load point required. A heading card for each set specifies 
the load point number and a scale factor to be applied to all forces in the 
set. The scale factor can be used to reduce the number of digits required 
to specify the forces, change the sign of the force and to conveniently 
scale the force without modification of the entire table. The heading card 
of each set is followed by cards specifying the forces in pounds at the pre- 
viously established times . Any number of tables may be read, with read-in 
of tables ending when a blank heading card is encountered. Any forces not 
specified will be treated as zero by the program. As each forcing function 
is read in, it is interpolated to a set of values corresponding to the time 
scale established on the time card and the resulting array is stored in a 
direct access data file for later use. 

Following the tables, a RESP card specifies which responses are to be 
processed and describes the filtering which should be applied to the ITF's. 

A PLOT card follows, which specifies whether or not the structural responses 
should be plotted. 
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For each response “to be evaluated, and for each load point, the required 
ITF is read from file 2 and unpacked. The ITF status array is read from file 
1. If no filtering is required, the zero adjustment value is obtained from 
the status array. If filtering is required, the ITF is passed to the FILT 
subroutine which removes specified low frequencies. The filtered or unfilte^ed 
ITF is now averaged over the time intervals specified on the time card and 
zero adjusted. The required forcing function is used from the direct access 
data file, and the Duhamsl integral is evaluated. Finally, scale factors 
for the ITF and the forcing function are applied. The result is the struc- 
tural response at one response point to one forcing function. The operation 
is repeated for each forcing function, with each response being printed and/or 
plotted, and accumulated in a summed response array which eventu all y gives the 
total response at the point to all the supplied forcing functions. The sunned 
response is printed and/or plotted. 

After all specified responses have been evaluated, a CYCLE card allows 
the program to recycle to different places in the input sequence to perform 
additional operations. The forcing function tables do not have to be repeated 
if they are unchanged. 

The filtering technique, described in Section 2, is applied in subroutine 
FILT. The sequence of operations closely follows that description. A feature 
of the filtering which. needs further explanation is the apparent elimination 
of frequencies from the high end of the spectrum as well as the low end. As 
an example, consider a case of 100 points at spacing of 10 ms, for a 1 second 
ITF. The fast Fourier transform fits this function with frequencies range 
from 0 cps to 99 cps in steps of 1 cps. In order to eliminate frequencies at 
or below 10 cps the first 11 and the last 10 coefficients would be set to 
zero. The last 10 are zeroed because the apparent fitting of 90-99 cps does 
not really evaluate this frequency content (there is only 1 data point per 
cycle). It does in fact produce a fit for the low frequencies 1-10 cps through 
a process known as aliasing. This is a characteristic of the fast Fourier 
transform technique, and is taken into account by setting to zero the higjh 
order coefficients as well as the low. 
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The ITF processing program performs seven functions in addition to pro- 
cessing Phase I tapes. Each operation is controlled hy a single input card, 
containing the function name and additional information when required. The 
functions, which are explained in the Program Description, are CLEAR, LOAD, 
DUMP, STATUS, LOADGQ, PLOT and END. Phase I processing is evoked by the 
function RUN and requires an additional data card. Tape unit logical numbers 
are under control of the input and should conform to the EXEC- 8 control cards 
for the run. 


Input Data 

Each operation is evoked by a single function card (2nd card required 
for RUN), with the following format. 

Function Card - format (A6, 1014) 

col 1-6 - function name, left adjusted 

col 7-10 - INN(l) 

i 

J appropriate integer data as described below 

col 43-46 - INN (10) 

CLEAR - no values of INN required. 

All data files are cleared. 

Usually performed before creating original data tape of a series, 
LOAD - INN(l) fortran unit number of data bank tape to be loaded. 

INN(2) number of first response on data bank tape (normally 1, 11, 
21, 31, etc.) This is used to check tape sequence number. 

DUMP - INN(l) fortran unit number of tape to which files will be dumped. 

INN(2) number of first response of the 10 to be dumped (l, 11, 21, 
etc.). This is used to determine the sequence number. 
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STATUS - INN(l) number of first response in a sequence, whose status is to 
he printed. 

INN(2) number of last response in sequence. It is assumed that 
appropriate responses are loaded in files. 

LOADGQ - INN(l) fortran unit number of tape to be loaded. This is a special 
tape for transferring data from the Grumman data bank to the 
MSC program. 

PLOT - INN(l), INN(2), INN(3), INN(U) 

Routine prints or plots ITF for responses INN(l) through INN(2) 
for loads (impulses) INN(3) through INN(U) 

INN(5) spacing for prints or plots in milliseconds 
INN(6) number of points to be printed or plotted 
INN(7) = 1 plot only 
= 2 print only 

* 3 plot and print 

INN(85 * 1 open plot routine for this set of plots 
= 2 close plot routine after these plots 

a 3 open and close plot routine 

This is used to initialized or terminate plot tapes as required 
by system. They are not opened or closed automatically by 
routine since several PLOT calls may be made on any one run. 

END - no values of INN required - terminates run. 

RUN - INN(l) fortran unit number of Phase I tape to be processed. 

INN(2) load (impulse) number 

INN(3) number of first response to be processed 

INN(U) number of last response (inclusive) to be processed 

INN(5) * 0 minimum number of records 

» 

! printed in output 
I 

= 5 maximum number printed 

RUN card must be followed by a single data card containing 
the full scale calibration values for the ten responses being 
processed, and the impulse track format (11E7.0) 
col 1-7 impulse calibration (lbs at full scale) 



col 8-lU response 1 cal. (units at full scale) 

i 

l 

i 

col 71-77 response 10 cal. 

Sample Data Setup 

Figure B-l shows a data sequence in which the primary operations 
are to process two Phase I tapes for the responses 1-10 for impulses 5 
and 6. 

The exact operations are: 

1 Load responses 1-10 from data bank tape mounted as unit 10. This 
has been created on a previous run. 

2 Print status of responses 1-10. 

3 Process Phase I tape on unit 11. It is an impulse 5 data tape - 
responses 1-10. Calibration data card follows. 

k Process Phase I tape on unit 12 - impulse 6, responses 1-10. 
Calibration data card follows. 

5 Print and plot the stored responses 1-10, impulses 5-6. 

6 Print status of responses 1-10. 

7 Dump files to unit 20 to create an updated data bank tape with 
impulses 5-6, responses 1-10 added or corrected. 

8 End run with normal termination. 

Sample Deck Setup 

A deck setup for the preceding data processing operations is shewn in 
Figure B-l. 

1 Fastrand assignment for temporary file 1, 2 tracks. 

2 Fastrand assignment for temporary file 2, 7 positions. 

3 Data bank tape on unit 10 XXX = tape label 
k Phase I tape on unit 11 XXX * tape label 

5 Position to file 2 

6 Phase I tape on unit 12 XXX = tape label 

7 Position to file 2 

8 New tape on unit 20, to be saved 

9 Fortran program ITF1 and subroutines with appropriate FOR cards 



10 Map cards 

11 XQT, data and FIN 

This setup does not include any control cards for plotting. 


27 


Output Data 

Primary printed output from the program consists of full listings of 
some records from the Phase I tapes, a summary of the average values of a 
member of records, listing of a tape status and printed ITF*s. 

Phase I Tape Records 

Figure B-2 shows a 3 page list of a record from a Phase I tape. The 
record covers 32 milliseconds, indicated on extreme left. The impulse track 
has 20 samples per millisecond and the responses have 5 samples/ms. Responses 
are always numbered 1 to 10 regardless of true response numbers associated 
with the tape. The numbers displayed are in counts which range from -2048 to 
+2047. Seme notable features of the record shown are: 

1. The impulse appears during the 26th ms of the impulse track and 
can be estimated roughly to .2 count-seconds 1.6 lb-seconds 
based on a calibration of approximately 12,500 lbs = 1600 counts. 

2. Response 11 is a dummy, 

3. Response 8 shows an oscillation of 16 ms period (60 cps), probably 
the result of electrical interference. It is removable by filtering 
but makes the estimated zero unreliable. Status output for this 
response shows a larger discrepancy between the zero derived from 
pre-impulse records and the zero obtained by ITF averaging than is 
normal . 

4. 43 numbers at the bottom are 32 - 1 ms averages for the impulse 
track, one 32 ms impulse average and 10 - 32 ms averages for the 
responses. 

A summary of record averages at the end of processing a tape shows the 
average value in counts of selected records for the impulse and 10 response 
tracks (Figure B-3). 
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Output from the STATUS operation presents the six items in the status 
file for the impulses and responses requested. The six items are: 

1* Conversion factor - engineering units /count. 

2. Zero adjustment - derived from region hetveen event pulse and 
impulse. 

3. Duration of stored response in milliseconds. 

4. Zero adjustment - derived from record average of response, omitting 
first 10 seconds (320 ms). 

5. Impulse scale factor - lbs /count. 

6. Impulse area in lb-seconds. 

Printed output from the PLOT operation is in zero adjusted counts as 
shown in Figure B-4 according to specifications on the PLOT function 
card. 

Plot 8 

Actual calls to plotter software have been omitted since plot routines 
under EXEC-8 were not wholly operational at time of programing . Further- 
more, final layout of plotted output cannot be properly determined until a 
specific application is undertaken. These plots are considered to be inter- 
mediate output, for checking purposes, and are best designed to match oscillo- 
graph playbacks of the original analog tapes. The arrays required for plotting 
are set up in the PL0T1 routine and correspond to the printed output. 
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The response program computes the structural response to a set of 
forcing functions, applied at one or more load points. The data required 
to produce one response consists of the time histories of the forcing func- 
tions and the set of ITF's relating the load points to the response point. 
Output from the program is printed or plotted structural response time 
histories . 


Input Data 

1. LOAD card FF, IU, IR1, IR2 (A**, 2 X, 31*0 
col 1-H FF. - ’LOAD' (literal) 

7-10 IU - tape unit on which data "bank tape containing ITF's is 
mounted 

11-1** IR1 - first ITF number on data bank tape 

15- 18 IR2 - last ITF number on tape (IR2 * IR1 + 9) 

2. TIME card FF, NPTS, NMS (A**, 2X, 21*0 
col 1-1* FF - •TIME* (literal) 

7-10 NPTS - number of points to be calculated in the resulting 
time histories 

ll-ll* NMS - spacing, in milliseconds of the response points 
NPTS = 25,0, NMS * 2 results in a time history of 500 ms - 250 points, 
every 2 ms . 

3. FORCE card FF, (MT(l), DT(l), I ■ 1, 6) 

(A5, 5X, 6(15, F5.0) ) 

This card establishes the time scale for the set of forcing functions 
which will follow. It assumes a set of values will be given at a fixed 
time increment, a second set at another time increment, and so forth, up 
to six sets if required (e.g., 50 points at 2 ms, 10 points at 10 ms, 

1 point at 100 ms). 

col 1-5 1 FORCE* (literal) 

11-15 MT(l) number of points in first set 

16 - 20 ET(l) increment in ms of first set 

21-25 MT(2) number of points in second set 
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26-J-30 DT(2) increment in ms of second set 

f 

76-70 DT{ 6 ) increment in ms of 6 th set 

If JfT(l) * 0, no forces are read 

4. Forcing Function Tables 

A set of force time histories for each forcing function must be 
supplied. The time scale for the forces is supplied by the FORCE 
card. The number of values required is equal to the sun of MT(l) 
on the force card. Forcing functions may be supplied for any or all of 
the load points for the problem. Those not supplied will be treated as 
zero. The table for each forcing function consists of one heading card 
followed by the required number of forces, 16 to a card. (If MT(l) 
from the FORCE card is zero, no tables are read.) 

a. Heading Card IL, SCL (12, F8.0) 

Col 1-2 IL load point number 

3-10 SCL scale factor to be applied to all forces in table. 

This may be used to change the sign of the force, or 
reduce the number of digits required in the tables (do 
not omit). 

b. Force cards (FEMP(I), 1 = 1, MTT) (l 6 F5.0) 

As many values of the force table as specified by the FORCE card, 

16 to a card in 16 F5.0 format. Force is in lbs. 

MTT 

c. End of tables card 

A blank card is required at the end of the tables to terminate the 
table read-in. 

5 . RESP card FF, IRA, IRB, IZFIL, FREQ (A 6 , 314, E12.0) 

This card specifies which responses are to be calculated and how they 
are to be filtered. 


^ MT(I) 
1=1 
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col 1-4 FF - ’RES?' (literal) 

7-10 IRA - number of first response in sequence to be calculated 
11-14 IRB - number of last response in sequence to be calculated 
15-10 IZFIL - 1 do not filter; use zero adjustment calculated from 

time before the impulse 

2 do not filter; use zero adjustment calculated from 
ITF record averages 

3 use Fourier Filter 

19-30 FREQ f vhen Fourier Filter is used, eliminate all frequencies 
below FREQ in cycles/sec* 

6. PLOT card FF, NPL (A4, 2X, 14) 
col 1-4 FF - ’PLOT’ (literal) 

7-10 NPL - 0 no plots 

1 plot total response 

2 plot individual response to each load as well as 
total response 

7. CYCLE card FF, NCY (A4, 2X, l4) 

After the PLOT card is read, the set of responses specified on the 
RESP card are evaluated and the results printed. The program may now 
be made to recycle to a point in the input stream to read data for a 
new set of cases with a minimum of data repetition, 
col 1-5 FF - ’CYCLE* (literal) 

10 NCY 1 return to beginning of data LOAD card 

2 return to TIME card 

3 return to FORCE card 

4 return to first force header 

5 : return to RESP card 

6 terminate program 

If the program recycles to the LOAD or TIME card, and the forcing 
functions are the same as for the previous run, setting MT(l) ■ 0 on 
the FORCE card will by-pass reading of the forcing functions and use 
the values from, the * previous case. 



Sample Deck and Data Setup 

Figures B-5, -6 shew a sample deck setup vith data. 

1. Fas tr and temporary data file assignments for Files 1, 2 and 3* 

2. Binary data hank tape containing ITF’s for response points 1 through 

10 . 

3. Fortran decks 

4. Map cards 

5. Load card specifying fortran unit 11, responses 1 through 10 (see 
item 2). 

6. Time card specifying the responses are to be calculated as 200 
points every 2 milliseconds, for a 4 00 ms duration, 

7* Force card specifying that the forcing function values are given 
at 2 ms intervals for 63 values (126 ms) with a final point at 
1126 ms. 

8. Forcing functions for 24 load points. Each set contains 64 
values (see item 7)- Scale factors for all 24 are unity, with 
slgxi changes to account for coordinate system definition. 

Forces are in lbs . 

9. Blank card to terminate table reading. 

10. RESP card specifying responses 1 through 4, to be run vith Fourier 
Filtering, eliminating frequencies below 5.2 cycles/sec. 

11. Plot card specifies emit plots. 

12. Cycle card specifies terminate run. 

Deck setup does not include control cards for plotting. 

Output 

Output from the program consists of tables of the forcing functions, 
interpolated to the time interval specified on the TIME card, and listings 
of both the responses to each individual forcing function and the summed 
response to the entire set of forcing functions. The units of the responses 
correspond to the units originally specified by the engineering units cali- 
bration values supplied in Program I when converting the Phase I tapes. 



33 


Plotting 

As in Program I, a plotting program is provided which does not contain 
actual calls to the plotting software. Arguments supplied to the skeleton 
routine are: IR, IL, MPT, NMS, Y. 

IR = response number 

IL * load number 

MPT = number of points provided 

NM5 * spacing of points in milliseconds 

Y * array contain MPT values of the response in appropriate 
engineering units - corresponds to the printed output. 
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5 . CONCLUSIONS & RECOMMENDATIONS 

Both the ITF program and Response program operate with a reasonable 
degree of efficiency in their present form (i.e., tape data base), provided 
the operations are sequenced to avoid excessive loading from and dumping to 
tape. This limitation could be removed by the use of a direct access storage 
device as the permanent data bank, as is the case for the Grumman system. 

While the program handles an unlimited number of response points, the 
number of load points is presently set to 36. This number can be easily 
increased by changing the DEFINE FILE statements, appropriate DIMENSION 
statements and the counter, NI, at the start of each program. 




FORfiS ITfltlTFl 

IMPULSE TRANSFER FUNCTION - ITF PROGRAM 1 
TO PROCESS PHASE 1 DIGITAL TAPES AND STORE ITF*S 

_ _MIN. PROGRAM _ _ 

COMMON I IN > 1 OUT *NR *NI • I LOAD 
IFILE DIMENSIONED TO 6*NT {MINIMUM OF 3000) 

DIM ENSION IFILEOOOO) _ — -• 

INTEGER FUN C # FUNCX 
DIMENSION FULL 7 ! 10 ) 

DIMENSION STAT (j6# 36 ) _ 

DIMENSION INN(iO) ♦FUNCUO) 

DIMENSION NTEST ( 8 ) 

DATA FUNC /'CLEAR '••LOAD * » * DUMP LvlRUN 1 * STATUS , _* 

" 1 • LOADGG* V ' PLOT ' r » * DUMMY 1 t 1 DUMMY > • * END * / 

900 FORMAT ( A6*10I4) 

901 FORMA T (//« FUNCT I ONXAflD^^UA^mJAj 

902 FORMAT ( // * CLEAR FUNCTION COMPLETE * ) 

903 FORMAT ! //* FILES LOADED FROM UNIT*#I4#» RESPONSES * I 3 ♦ 

904 FORMAT!//* FUNCTION NOT I DENTIFIED * ) 

90 5 _ FORMA T ( //* FILES DUMPED TO UNIT 1 *14* 1 RESPONSES ‘t 1 3 * » « 

906 FORMAT ( HE? .0 ) 

907 FOR MAT ( 1 H X / [ _R£,S,PQNSE- NO. * *, 14// j 

908 FORMAT !I4*6F12#4) 

909 FORMAT!//* ERROR IN INTEGER FIELDS* *414) 

910 FO RMAT! 1 HI / * EXECUTE PR OGRAM I - ITF PROCESSING*,.). 

911 FORMAT!//* ERROR* TAPE MOUNTED IS IS FOR SET* #14) 

912 FORMAT!///* TERMINATE RUN*) 

91 3 FORMAT ! 2 16 F 10^6 ) . . 

' 914 FORMAT ! 3000 1 6 ) 

915 FORMAT!//* PLOTTED I MPULSES * # 14 ♦ *- * * 1 3 * * RESPONSES ** 14 » 

JCV!J)«U-1)/NR+1 — 

CALL MFC (NTEST) 

SET I/O DEVICES 

I IN* 5 _ — 

rour-6 

SET FILE SIZE 

* NUMBER OF RECORDS IN FILES-NT _ _ 

: prun nr records of 6*Ni words 

: FILE 2 10#NR*NI RECORDS OF 200 WORDS 

DEFJNE FILE 1J 10*216* Uj_IUl ) _ „ 

"DEFINE FILE 2 t 3600*200 *U#fU2") 

• NR IS NO. OF RESPONSES PER TAPE 

: N I IS NO., OF IMP ULS ES IN. . SY STEM 

NR* 10 
N I * 3 6 

NT°N R«NI _ _ 

NTI*6#NI 

NTT*6*NT 

nttt*jo#nt - 

"12.0 

WRITE! IOUT*910) 

: READ FUNCTION CARD_ . ; — . 

X READ < 1 IN *900) FUNCX * INN 

WRITE! IOUT #901 ) FUNCX * INN 

: CHECK FUNCTION _ ... . 

DO 4 I FUN C * 1 # 1 0 


" THR00010 

THR00020 
THR00030_ 
THR00090 
THROOIOO 
THROOUQ 

THR00120 

THR00130 

THR00140 

THR001 50 . 
'' THR00160 

THR00170 

THROO lQO 

THR00190 
-* * 1 3 ) THR00200 

THRQ-O-ZjLQ- 

* *13) THR00220 

THR00230 

IHR.00240__ 

THRQQ25Q 

THR00260 

THR 0Q270. 

THR00280 

THR00290 

thro 03_oo 

THR 00310 
-* * 13 ) THR00320 

THR00330... 

THR00340 

THR 00350 

THRbb360 ' 

THR00370 

THR00040 

THRQQQ5Q 

THR00060 
THR00070 
THRbbbeb 
THR00380 
XHR00390 
THR00400 
THR00410 
THR0P42_0_„ 

THR00430 
_TH R00440 
THR00450 " 
THR00460 
THR00470 

THR00480 

THR00490 

111RQ0500 

THRobsio 


Figure A-l 

Program I - Fortran Listings 
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1F<FUNCX«EQ»FUNCUFUNC) ) GO TO 6 . . . - 

4 CONTINUE 
ILLEGAL FUNCTION 

WRITE ( IOUT #904) . 

THR00520,. 
THR00530 
THR00540 
THR005S0 

GO TO 1000 
PERFORM FUNCTION 

6 GO TO ( 100*200 #300 #400 #500 #600 #700 #800 *950»1000)#I FUNC ...... 

THR00560 

THR00570 

THR005§Q_~ 

CLEAR BOTH FILES 
100 CONTINUE 

DO 7 J*1 # 10 

THR00590 

THR00600 

7 WR1TE(1'J)<IZ*I»1»NTI) 

DO 102 J-l.NTTT 

102WRITE<2'J)(IZ«I«1*200) .... ... ... .. 


WR ITE < IOUT #902 ) 

GO TO 1 

LOAD BOTH FILES FROM TAPE 

THRQ0640 
THR00650 
THR00660 

200 CONTINUE 
IN1»INN( 1 ) 
IN2*INN< 2 ) 

THR00670 

THR00680 

THR00690 

IN3-IN2+9 

CHECK TAPE SEQUENCE NUMBER 
I CODE" ICV ( I N2 ) 

THR0Q700 

THR00710 

THR00720 

REWIND INI 

READ ( INI ) (IC0D#I*1*100) 

1F( t CODE •E0» I COD ) GO TO 201 . 

THR00730 

THR00750 

WRITE! I0UT»91 1 ) I COD 
REWIND INI 

GO TO 1000 

THR00760 

THR00770 

THR00780 

['201 READ tINl)(lFILEli)*I*l*NTTI 
DO 203 J«l»10 

11«NTI*( J-lJ+1 .... .. - ... 

I2-I1+NT 1-1 

203 WRtTEtl'J) (IFILEUI *l«Il*I2) 

DO 202 J*liNTTT *10 

THR00810 

READ UNI) UFILEU ) #1*1 #2000) 

DO 202 JJ«1#10 

LL«200#JJ ... ... 

THR00820 

L*LL-199 

JJJ«JJ+J“1 

202 WRITE ( 2* JJJ) ( IFILE(K) »K*L#LL) _ . 

, REWIND INI 

WRI TE ( IOUT #903 ) IN1»IN2#IN3 

' GO TO 1 ' 

THR00840 

THR00850 

THR00860 

DUMP’ BOTH FILES TO TAPE 
300 CONTINUE 
INI* INN ( 1 ) 

THR00870 

THR00880 

THR00890 

!N2*1NN(2) 

IN3-IN2+9 

rewind ini ...... ... .... . 

THR00900 

THR00910 

THR00920 

' WRITE TAPE SEQUENCE NUMBER 
ICOD“ICV< IN2) 

! WRITEUNI) UC0D#I*1#100) 

THR00930 

THR00940 

THR00950 

\ - “ DO 3^1 T*1 #10 ’ 

l»i*nti 

\ K.L-NTI+l ; 

r$ o'l' R EADlin ) TiF I LETJ) #~J«K#L) 
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WRITE! INI MIFILEUl.lFliNTT) ... . — 

DO '302“ J*1 *NTTT • 10 
00 303 JJ»1 1 10 

LL*200#JJ . . 

THR00980 

L*LL-199 

JJJ*JJ+J-1 

303 READ ( 2 * JJJ) ( IFILEt-K ) *K“l*LL) 


302 WRITE (INI) UFILEU) ♦ 1*1 *2000) 

REWIND INI 

WRITE(IOUT»905) INl*IN2*IN3 

THR01000 
THR01010 
THRO 102Q__ 

GO TO 1 

PROCESS A PHASE 1 TAPE AND STORE IN FILES 
400 CONTINUE ... 

THRO 1030 
THR01040 
THR01050__ 

ITAPE»INN(1I 
I LOAD* INN ( 2 ) 
IRESP*INNI3) 

THRO 1060 
THR01070 
THR01080 

IRESPP*INN(4) 

1TEST*INN (5 ) 

IRX-IRESP-! t JRESP-1 ) /NR)»NR 

THR01090 

THR01100 

THR01110 

IRXX*IRX+IRE$PP-IRESP 

IF(INN(2).LT*1»0R*INN(2)*GT*NI) GO TO 440 
READ ( 1 IN* 906 ) FULL! »FULL 

THR01120 
THRO 1130 
THRQ1140 

CALL STORED! ITAPE * I RX* IRXX t ITEST *FULLI *FULL*NG) 
GO TO X 

440 WRITE(IOUT*909> ( INN! I) ♦ 1*1 *4 » 

THR01150 

THR01160 

THR01170 

GO TO 1 

PRINT STATUS OF TAPE 

500 CONTINUE - - 

THRQ1180 
THRQ1190 
THRO 1200 

INI* INN! 1 ) 

IN2-INNI2) 

DO S10 IR*INltIN2 

THR01210 

THR01220 

THR01230 

IRX*IR-( ( IR-i)/NR)*NR 
WRITE! TOUT. 907) IR 

READU'IRX) ( (STAT! I «J) *1*1 *6) *J*1*NI ) _ 

THR01240 

THR01250 

DO 506 IMP*T*NI 

506 WRITE( IOUT *908 ) IMP* (STAT ( I • IMP) *1*1*6 ) 
310 CONTINUE 

THRO 1270 
THR01280 
THRO 1290 

00 TO 1 

LOAD FILES FROM GRUMMAN BCD TAPE 
600 CONTINUE 

THRO 1300 
THR01310 
THR01320 

INI* INN! 1 ) 

CALL MRWND(INl) 

DO 610 1*1*10 __ . _ , 

THRO 1330 

' j A L L MREADUNl.lVlFILE.360.NTEST.MRD) 

611 CONTINUE , 

IF! NTEST ( 6 > •EQ»1> GO TO 611 ■ 

' " "CALL CRONCHTTFILE .360) 

WRITE! 6*930) (I FILE! J).J«1*360) 

o BflRMAT!lHl/10!lX*012)l ■_... 

DO ~ 613 K“1 #36 
KK*10#K-9 

DECODE (612#IFILE(KK) ) t STAT ( J*K ) * J*1 * 6 ) 


612 'FORMAT! 6F10. 6 ) 

613 CONTINUE 

DO 710 1FIX*1 *36 — 


710 STAT (3 *1F IX) *3000. 


Figure A-l (Cont f d) 


Program I - Fortran Listings 


39 





IJl<KS.LT*KSX+20) GO TO 10 

OVER 6 SECONDS OF DC CALS 
KFAIL*2 

fiSLJLSLML ...... 

34 IRECY-IREC 

k$y*k$-i 

I.M-3 - ... 

GO TO 10 


FOJJOQ52Q 
FOU00530 
FOUQ0540 
FQ.UQ0550 
FOU00560 
FOU00570 
. FOUPOSSO 
FOUQ0590 
FQU00600 


C IN. ZERO ARM, 

ffi SKIP 50 RECORDS 
■ 40 IF I KS-KSY*NE* 5 ) GO TO 10 




SEARCH 150 RECORDS FOR EVENT PULSE 
KAVG*I AVG ( l'»KS) + 120 


K$*KS+1 

IREC-IREC+1 


NR£CXKSJ..»1R£C 

NG»3 

CALL READ7 ( ITAPE# IREC1 » IREC2 » I AVGtKStNG ) 


_ __LF1NG « EQ.1 )„ GP„T 0_9 9_ 

IF( IAVG( 1 iKS ) *GT«KAVG) GO TO 54 
52 CONTINUE 

FAIL TO FIND EVENT PULSE. 

KFAIL-3 
GO TO 99 


F0yp06l0.__ 

FOUQ0620 

FOU00630 

F_Qy00640L_ 

FOU00650 

FOU00660 

F0UQQ6 70 

FOU006Q0 

FOU00690 

.Fjouoj37iia..: 

FOU00710 

FOU00720 

FQU0 0730 .... 

FOU00740 

FOU00750 

F OUO Q76Q 

FOU0Q77O 

FOUOO70O 



FOUND EVENT PULSE RECORD 
54 CONTINUE 

IRECZ«IR£C . 

KS2-KS 

READ__6 RECORDS AND SAVE LAST I N ENTIRET Y IN IDISK ARRAY 
‘DO 60 1*1 #6 
IREC-IREC+1 


FQUQ0790 
FOU00800 
FOU00810 
FQUQ0 820 
FOU00830 
FOU00840 
FOU00850 
FOU 00860 
FOU00870 


NG* 3 . „ _ 

K$*KS+1 

CALL READ7UTAPE* IRECl # IREC2 # I AVG# KS#NG > 


l 


IF(NG.EQ .l) GO TO 99 

NREC (KS ) * IREC 

60 CONTINUE 

IFJJ TEST «GT«0 ) _ 

1CALL RITRECC IREC* IREC2 • IREC2 < 641 ) ) 
DO 61 1*1 *640 

61 IREC22U )*IREC2(I) 

DO 62 J*lTiO 

DO 62 1*1 #160 

62 IDISK! 1 4200 # J ) * JJREC t I ♦ J ) _ __ 

READ AND SAVE ANOTHER RECORD 

JREC-IREC+1 

K$*KS+1 


FOU0Q8 80 

FOU00890 

FOU00900 

F0U QQ91Q 

FOUQ0920 

FOU00930 

FO U 00 9 40^ 

FOU00950 

FOU00960 

FOU009J70 

FOU00980 

FOUQ0990 

FOUQ IOOQ 

fouoioIo 

FOU0102 0 
FOU01030 
FOU01040 
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KSF-KS 

NG-3 

CALL READ7 ( I TAPE* IRECl * IREC2 » I AVG#KS*NG) 

IF (NG»EQt 1 ) GO TO 99 _ 

IF ( I TEST • GT*0 ) 

KALL RITRECC IREC* IREC2* 1REC2I641) ) 

NREClKSJfjREC _ _ _ 

00 63 1*1 *640 

63 IREC22U+640)«IREC2(I > 

DO 64 J*l*10 

64 I * 1*16 0 

64 IDISKU+360*J)«JREC(I#J) 

“SEARCH 2 S^E D^ ^CORDS~ FOR TMP U L SE 

lAVGX*<IAVG(l*KS-2>+IAVGU#K5-3) + IAVGU*KS~4n/3. 

DO 65 1*3*12 50 : 

IF( IREC22 ( I ) *GT* I AVGX+200* «AND« IREC22 ( 1-1 ) • GT . I AVGX+200* • AND. 
1 IREC22U-2UGT.IAVGX+200.) GO TO 66 
63 CONTINUE 

IMPULSE NOT FOUND 

KFAIL»4 _ _ ; 

GO TO 99 

FOUND IMPULSE I * THIRD ITEM OVER 200 

COMPUTE AREA OF IMPULSE 

GET LOCAL ZERO „ F^_1MP ULS_L AREA CALC ULATION 

66 J*I-2 
XI A$H*0« 

K*J-43 _ _ _ _ 

IFU.Lt.l) GO TO 99 
KK-J-6 

DO 67 L*K *XK 


67 XI ASH* XI ASH+IREC22 ( L ) 

XtASH-XIASH/40* 

I A* X I ASH+SIGN ( •5* XI ASH ) 

FIND FIRST AND LAST POINT 50 OVER AVERAGE 
66 I F < IREC22 ( J“1 ) «LT • IA+50 ) GO TO 70 


Jm 

GO TO 68 

__7^J^iIREC22U + lu>*iJ*IA+J.01 GO TOl_H 
I-I + l 
GO TO 70 

7 1 JSTRT *J . _ 

CALCULATE AREA OF END PIECES 

AREA^000025#< (IREC22< JJ-X I ASH ) #*2/ U REC22 < J >-I REC22 < J-l ) ) + 

‘ T'lIREClKTj-XIASH 1**271 1 R EC 2"21T ) - 1 RECT2 ( I + 1 f )T 
IPTS-I-J+1 

IFU IPTS/2)#2-IPTSiNE.0> GO TO 74 
TOUffTTW" ODDNuMBER OF 'POWS" FOR SIMPSON RULE “ 
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+ AREA 


J-J+X 

AREA- • 000025* ( I REC22 1 J > + IREC22 < J-i ) -2 **X I ASH ) 
PERFORM SIMPSON INTEGRATION 

74 tSUM«IREC22( JH1REC22U ) „ 

75 JX-J+1 
IX-I-X 

iFUitLT.Jl) GO TO 99 __ _ 

00 76 K-JX •11*2 “ 

76 ISUM«I$UM+4*IREC22 HO 
JX-JX+X 

rx^rw 

IF(n.LT.JX) GO TO 78 
DO 77 K-JX * I X * 2 

77 ISUM-ISUM42#IREC22( K) 

78 AREA-AREA*. 00005/3»*ISUM-.00005#<i-J> *XlASH 


COMPUTE CALIBRATION VALUES - ' 

IFUSY-KSX.GT.4) GO TO 79 

NOT ENOUGH DC CAL RECORDS 

KFAIL-5 
GO TO 99 

79 FULLI-3.» FU LL I / (I AVG U tKSX + 2 ) + IAVG ( X »K5X+3 H-IAVG 

X -IAVG<X*XSY+X) ~I AVG { 1 *KSY+2 ) -I AVGU »KSY+3 ) ) 


( X • KSX+4 ) 


AREA IN PO UND-S E CONDS .. 

AREA Z-AREA*FULLI 
DO 42 ! -2 iXX 
J-t-X 

42 FULL t J ) -3 .*FULLC J ) 7 IIAVG( I *KSY-3 )+IAVG< I tKSY-4)+I AVGU *KSY-5 ) 
X -I AVG< I »KSY+X ) -I AVG 1 I #KSY+2 > •*! AVG( I iKSY+3 ) ) /AREA I 


IT IS NOW NECESSARY TO FILL THE - ZD I SK ARRAY TEN TIMES AND 
IT INTO FILE 2 TO COMPLETE STORAGE OF RESPONSES 


WRITE 


STORE 


TSWnWTAVEb IN I DISK ARRAY' 
LAST KNUM MILLISECONDS OF THESE 


VALUES 


K NUM *Tl 2 9 6 * J S f R f ) 7 2 0 
K-3X6-5*KNUM 
DO 228 I -X tKNUM 
K-K*5" 

DO 228 J-X*XO 

AVERAGE 5 P OINTS TO GET X MILLISECON D DATA 

228 IDISMI t J)-"< ID Is‘k (K + 200*J)-M DISK (K + 201 *J) + I DISK (K+202» 
X IDlSK<K+203*J>+IDISKtK+204»J>+2>/5 

NSET* X _ _ ___ ... 

I -KNUM 

K-X56 

DUR-60 00* 

GO TO 211 
209 DO 2X0 J-XtXO 

2X0 1 DISK 1 1 • J ) ■ ( JREC( K* J)+JRECUtX_t^)+JRECUt2|JKURECjK+3 
1 “ JR EC ( K+4 1 J ) +2 ) /5 


• J) + 


_ FOUO1580 
FOU01590 
FOUOX600 
F0UQ161Q 
FOUOX620 
F0U0X630 
F0U0X 64Q 
F0U0X650 
F0U0X660 
F0U01670 
F0U01680 
F0U01690 
FOUOX700 
F0U01710 
F0U01720 
F0U0X 73Q 
FOU01740 
F0U0X750 
F0U0 1760 
F0U01770 
F0U0X780 
FQU0 X790 
FOUOX800 
FOUOX8XO 
FQU0X820 
FOU0X830 
FOU0X840 

F QU0X8 50 

FOU0X860 
F0U0X870 
__FOjJ01880 
FOU0X890 
FOU0X900 
F0U 0X9XQ 
FOUOX920 
FOUOX930 
FOU0X94O 
F0U0X950 
FOU0X960 
FOUOX970 
FOUOX986 
FOU0X990 
_F0U02000 
FOU020X0 
FOU02020 
FOUQ 2Q3Q 
FQU02040 
FOUQ205Q 
_ _FOUO_2060 
FOU02070 
FOUG208Q 
F0U Q209Q 
F0U02X00 
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_H1_IPU«NE«156> GO TO 212_„ _ 

K*-4 ~ 

KS-KS+1 

IREC*! REC+1 

NG*2 

CALL READ7 (ITAPE* IREC1 * IREC2 *1 AVG*KS*NG) 

MKlKUfJRlC. 

!F<NG#NE*1) GO TO 212 
EARLY END TO DATA 
DUR*600«<NS£T-1)+ I 

n«i+i " ~ ” ' 

IF ( 1 1 #GT *600 ) GO TO 220 

K$M*KS-1 

DO 213 1*11*600 
DO 213 J*l#10 

FILL OUT RECORD WITH AVERAGE OF LAST RECORD 
213 10ISK(1#J)*I AVG ( J+l *KSM ) 

NSET*10 
GO TO 220 


FPU02110 

FOUQ2120 

FOU02130 

FOy0214Q 

FOU02150 

FOU02160 

FQUQ21T0 

FOU02180 

FOU02190 

FOU02200 


F0U02210 

FOUQ2220 

FOU02230 


FOU02240 

FOU0225O 

FOU02260 


4 


212 K*K+5 

IFU.EQ.600) GO TO 220 
1*1 + 1 

GO TO "209 ~ “ 

PACK 600 WORDS INTO 200 
220 L-l 

Ml*4096*4096 
DO 222 LL*1 *200 


FOU02270 

FOU02280 

FOU02290 


FOU02300 

FOU02310 

FOU02320 


FOU02330 

FOU02340 

FQ U02350 

FOU02360 

FOU0237O 

FOU02380 


DO 221 J*1*10 

221 !DISK(LL»J)*UDISK(L»J> +2048 > *Ml+( 1DISK ( L+l * J 1+2048 >*4096 + 

_1 < ID ISK < L+2» J) +2048 ) 

222 L*L+3 
JJ*0 

C STORE IN FILE 2 


DO 242 J* IRESP • I RESPP 
JJ*JJ+1 

IFILE*360»{J-1)+10»( ILOAD-l)+NSET 


FOU02390 

FOU02400 

FOUQ241Q 

FOU02420 

FOU02430 

FOUQ2440 


242 WRITE* 2 » I FILE MID ISK *I*JJ) #1*1*200) 
1*1 

NSET *NSET+1 

1 f ( n srTveo. rrr“Grro~T43 

GO TO 209 

349 WRITE* I0UT*107) AREA I *DUR 


107 FORMAT *777 • IMPULSE AREA I N TB^EC^TF 1 0 *6 / 

1 » DURATION OF RESPONSE IN SECONDS* 1 *F7 *0 \ 


FOU02450 

FOU02460 

FOUQ247Q 

FOU02480 

FOU02490 

FOU025OO 

T5U025IO" 

FOU0252O 

F0U02530 


CALCULATE ZERO FROM REGION AFTER EVENT PULSE 
DO 350 1*1*10 

l* i+i _ ; 

350 ZEROTl ) *T| AVGTL • KSZ+21 +TAVG t L * KSZ+3 ) + 1 AVG ( LVKSZ+4 > + 1 AVG *~L »KSZ+5TT 
1 /4. 

CAraJtATE“ALTERNATf~"ZE RO~FROR~ AVERAGrOF ENTn?E"TTF “EXCEPT 


FOU0254Q 

FOU02550 

FOU02560 


F0U02570 

FOU02580 

F0UQ 2590 

FOU02600 

FOU02610 

FOU02620 

75UOXSJO~" 
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JtA.ST.JRE CORO AND FIRST 320 MILLISECONDS 
K$FF«KSF+iO 
KSL«KS-1 


DO 33 4 I»l>10 

2ER0X ( I ) «0 • 

DO 332 K*KSFF tKSL 

3l5I_2E.RQX ( IJ-ZEROX LH + I AVG( I+1»K) 

^ 354 ZEROX ( I ) *ZEROX ( I ) / < KSL-KSFF+1 ) 

C UPDATE ST ATUS FILE„ 

i»o ‘ 


FOU02640 

F0UQ2650 

FOU02660 

FOU02680 

FOU02690 

FOU027PP 

FOU02710 
_FOJJ02 72_Q 
FduQ27~30 



DO 360 J»IRESP* IRESPP 
1*1 + 1 

FOU02740 

FOU02750 


READ < 1 * J ) STATUS 
STATUS ( 1 » ILOAD ) -FULL (I! 
STATUS (2 tILOAD) -ZERO (I) 

FOU02760 

FOU02770 

FOU02780 


STATUSO* ILOAD) "DUR 
STATUS (4 » ILOAD) »ZEROX( I ) 
STATUS 1 5 1 ILOAD) bFULLI 

FOU02790 

FOU02800 

360 
99 

STATUS 16* ILOAD) *AREAI 
WRITE! 1 ' J ) STATUS 
CONTINUE 

FOU02820 

FOU02830 

FOU02840 

101 

CALL MRWND C I TAPE ) 

WRITE <!OUT*101MJtJ*l#9) 

FORMAT < ///• SUMMARY OF RECORD AVERAGES*// 1 RECc IMP »t 



1 9( * R * » 1 1 ) # * RIO*//) 

Jl*KSX-4 
J2-KSY+4 

FOU02880 

FOU02890 


WRITE! IOUT. 100) (NRECIJ) »( IAVG! I »J) . I«1 .11 ) • J-Jl *J2 ) 

J1"KS2“2 

J2*KSZ+6 

F0U02 


WRITE! IOUT. 108) (NREC!J) .( IAVG! I *J) « 1*1 • 11 ) .J«J1.J2> 

Jl-KS-4 

J2-KS 

FOU02940 

F0U02950 

CD 

O 

H 

i 

WRITE(IOUT»108) { NREC ( J ) t ( IAVG U t J) *1*1*11) • J«J1 * J2 > 
FORMAT <1216) 

RETURN 

FOU0297O 

FOU02980 


END 

FOU02990 
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_f£38.»IS RITREC.RITREC 

SUBROUTINE RI TREC < JB . 18 . IC) 
COMMON I I N* IOUT »NR »NI *LOAD* I RESP 
DIMENSION ICSUMUll .IBSUM02) 

FI V0002 
FIV0005 

100 

DIMENSION IC(160.11).IB( 2400 ) 

FORMAT! '1RECORD NO. < .I4.10X. 'LOAD NO. ' .I3.8X. 'RESP. NO. '.13 
1 //40X.' IMPULSE TRACK' .60X. 'RESPONSE NO. 1'/) 

FIV0006 

FIV0007 

FIV0008 

101 

102 

103 

FORMAT (I3#2X#20I5»2X«5I5) 

FORMATtlHO # 1X#5 ( 5X ♦ 1 RESPONSE N0*'fI3*6Xn 
F0RMAT(I3#5I5#4(1X#5I5) ) 

FIV0009 

FIV0010 

FIV0011 

104 

105 

FORMAT ( 1H1 ♦. 1X*5 1 5X t 1 RESPONSE NO.*#I3#6XM 
FORMAT (lXtl6I5/lX#l6!5/I6#110»l0I5) 

DO 9 1*1 *11 

FIV0012 

FIV0013 

FIV0014 

8 

tCSUM( I ) *0 
DO 8 J*1 » 160 

ICSUMt I) -ICSUM (I)+IC(J# I ) 

FIV0015 

FIV0016 

FIV0017 

9 

ICSUMI I ) -ICSUM ( I ) /160.+.5 
WRITE! IOUT. 100) JB.LOAD. IRESP 
DO 4 1-1.32 

FIV001Q 

FIV0019 

FIV0020 


IB$UM( I ) *0 
J1*20#I-19 
J2*J1+19 

FIV0021 

FIV0022 

FIV0023 

12 

“D0~ir"X*~ J17 31 ’ ~ ~ 

IBSUM( n -I0SUMU )+lB( J) 

IBSUM( I)*IBSUM(I )/20#+#5 

FIV0024 

FIV0025 

FIV0026 

4 

J4- J3+4 

WRITE! IOUT. 101 )I »( IB(J) .J-J1.J2 ) .( IC< J.l I .J«J3. J4l 

FIV0027* 
F I V0Q28‘ 
FIV0029 


WRIT E< IO'UT • VD2T - Tx * k^76 ) 

DO 3 1*1*16 

Jl«5*I*4 

FIV0030 
FIV0031- 
FI VQQ32 1 

5 

J2-J1+4 

WRITE (IOUT #103) I * ( UC ( J #K ) • J*J1 # J2 ) #K*2 *6 ) 
WRITE (IOUT #104) (K#K*2#6) 

FIV0033- 

FIV0034* 

FIV0035 


DO 1 F*lf #32 ' 

J l*5#I-4 
J2-J1+4 

F I V0036 r 
FIV0037 f 
FIV0038- 

7 

WRITE! IOUT. 103 ) I . ( ( IC( J.K) • J* J1 . J2 ) .K-2 .6 ) 
WRITE! IOUT. 102) (K. K-7.11) 

DO 6 I >1 .32 

FIV0039 

FIV0040' 

FIV0041f 

6 

Jl*5#I-4 — — 

J2-J1+4 

WRITE! IOUT #103 ) I • ( ( IC t J »K > #J» J1 # J2 ) #K»7 #11 ) 

FIV0042 

FIV0Q43 

FIV0044 

LL 

ICSUMt 11 )*0 
DO 11 I *1 #32 

ICSUMt 11 )*IC$UM( 11) +1BSUM( I ) 

FIV0045 
FIV0046 
F I VGQ47‘ 


ICSUM! 11 ) * ICSUMt 1 D/32. + .5 

WRITE! IOUT. 10 3) I BSUM. ICSUM ( 11 ) • ( ICSUM (J ) » J«1 .10 ) 
RETURN 

FIV0048 

FIV0049 

FIV0050 


END 

FIV0051' 
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FIV00950 

FIV00960 

FIV 0 Q 97 Q 

FI V00980 

FIV00990 

FIVOiOOO 

FIVOIOIO 

F I VO 1020 

F I VQ 1 030, 

F I VO 1040 

FIV01050 

F I VQ1 06Q_ 

FIV01070 

F I VO 1080 

FIV01090 



K*641+(J-2 1*160 

L»K+1»9 

DO 4 I*K*L 

FIVOUOO 

FIV01U0 

FIVG1120 

4 

ISUM-ISUM+IRECU) 

FI VOX 130 

6 

IAVG(J)»ISUM/160.+.5 

FIV01140 


RETURN 

FIV01150 


END 

FIV0U60 



♦FOR# I S UNPACK #UNPACK 

SUBROUTINE UNPACK ( 1 N1 ♦ IREC2 # ISTAT ) 

DIMENSION IREC2 ( 2400 ) »N ( 8 ) #DUM(812) 

CAUL MFC(N) 

CALL MREADI I N 1 » 1 # DUM *8 12 *N »M ) 

3 C ONTINU E 

I FTNT6T7EQV1T GO TO r 
ISTAT ■ N<6) 

I R EC 2(1) *FLD ( 12 * 12 #DUM ( 1 2 ) > 

IR EC2( 2 ) ■ FLD (24# 1 2 # DUM (12)) 

IREC2 ( 2400 )*FLD( 0 #12#DUM(812 ) ) 

DO 2 K*13 #811 

IREC2( J)"FLD(0*12#DUM(K) ) 
IR£C2 ( J+l)ttFLD(1 2 »12 #DUM ( K ) ) 

2 I R EC 2 ( J+2 ) «F L D ( 24 #12#DUM (KM 
RETURN 
END 



*fPR#J$ AVG #AVG _ 

SUBROUTINE AVG UREC # IAVG #NG ) 
DIMENSION IAVGU1) 

.DIMEN SION I RE Cl 2400 ) 

t SUM-0 
DO 1 !-l#ll 

1 lAVG m«0 

C AVERAGE IMPULSE NG«1 

IF ( NG*EQ* 2 ) GO TO 3 

9 DO 2 I -1 #64 0 

2 ISUM* ISUM+IREC ( 1 > 

ISUM* ISUM/640 • + * 5 

IAVG i 1 ) -I SUM 

IFING*EQ.1) RETURN 

3 DO 6 J*2#ll 
I SUM-0 


RESPONSES NG-2 


r 



«FORtI$ PLOTltPlOU . 

S UBROUTINE PLOT 1 .... 

COMMON I 1N»I0UT#NR#NI 
DIMENSION STAT ( 6* 36 ) 

DIMENSION X ( 301 ) tV(301 ) » IFILEU ) »INN (_8J. 
IR1* INN ( 1 ) 

IR2* INN( 2 ) 

3J 

IL2*INN( A) 

MS*I NN ( 5 1 

NPT*INN ( 6 ) _ 


IPRMNNI7) 
IOPEN*INN ( 8 ) 

I F < IOPEN»EQ 1 1 >.OJ5,». 
XL* 10 • 


OPEN P LOT ROUTINE 


\ YL*2 • 

l X( ll»0i 

Yhl»0i 

, WRITEUOUTtlOO) 

100 FORMAT ( ///♦ PLOT PAC KAGE EXECUTIN G 1 ) 

h 1 NPTS-MS/NPT 

IF ( NPTS§LE«300 ) GO TO 2 

[ NP fNPT + 1 _ 

r GO TO 1 

2 DX-XL/NPTS 

[ DO 3 I-1»NPTS 

3 X<I+1)*DX*I 
MS*NPT*NPTS 

00 200 IL* IL1 « 1L2 

r DO 200 I R* IR1 # IR2 

' IRX-IR-M IR-n/NRMNR 

READ ( 1 1 1_RX ) t (STAT ( I ♦ J) t I*l»6 ) < J *1»NI ) 

[ C0N«STAT ( 1 1 1 L ) 

I 2ERO-STAT(2#IL) 

IF»{ 1RX-1 )»10*NI + 10* ( IL~ 1 ) 

.... ~ 00 6 J*l»10 

LL*200*J 

1 L-LL-199 : 

tFU3#L-3) .Gf .MS) GO TO 7 
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IFF* IF+J - 

6 READ(2*IFF) i IFILEIK) iK*L»LL) 

7 NPTSS-NPTS+l 

DO 40 1»1 #NPI5 

IP-I+1 
YUP J-0. 

DO 10 J» 1 #NPT _ 

K«(t-1)*NPT+J 
KW-U-U/3 + 1 

<P_M A0§ UC-1K W-l)»3 -l)«12 

10 Y 1 1 P ) ■ Y ( IP)+FLDIKP#12#IFILEIKW) ) 

40 YUPJ-YUPJ/NPT-2048.-ZERO 

_______ I F < I PR« EQ * 2 ) GO TO. 60. 

C PLOT X VS Y NPTSS VALUES _ 

C X IS IN INCHES RANGE 0 TO 10 INCHES IS TIME FROM 0 TO MS MILLISECONDS 

C __Y J_S IN COUN TS F RQM ... - 2048 TO +20. 4 8. — 

C CON IS MEASUREMENT UNITS PER PUUND SECONDS PER COUNT 
IF ( I PR« EO »1 ) GO TO 200 

60 WRITE<IOUT#902> I R » I L»NP TS » NP T »NPT _ 

902 FORMAT ( lHl/ 1 RESPONSE 1 *# 13 • * LOAD 1 #13# 10X # 14# ' POINTS# STARTING'# 

1 » AT ' # 1 3 ♦ * MS# SPACED EVERY'#I3#' MS IN COUNTS# ZERO ADJUSTED* ) 

WRI TE ( 1 0_U T #903 >J YJ I ) # I «2 »NPTSS) - 

903 FORMAT! ////( lH0#'ibFl0.1 > ) 

200 CONTINUE 

C IFdOPENiGTtl ) CALL CLOSE PLOT ROUTINE 

RETURN 

END 
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* FOR* I_$„R ESP • R £ SP _ 

DEFINE FILE 1 ( 10* 216 *U * IU1 ) 

DEFINE FILE 2 l 3600* 200 *U * IU2 ) 

DEFINE FILE 3 ( 36 » 12 00 »U * I U3 ) 

COMMON I IN* IOUT*NR*NI * ! LOAD 
INTEGER F #FF 

DIMENSION IF I L E ( 6000 ) » I TF ( 1200 ) *Y,lI,Fi 12,001, — — - 

‘ DIMENSION F(6)*RESPI1201> #SC<36) *MT<6) *DT(6) *TEMP< 1201) 
DIMENSION FEMP { 1 200 ) *JTEMP( 1200) *F 1 ( 1 200 ) *FOR ( 1200 > 

DIMENSION SUM (1201) *STAT ( 6 * 36 ) * I F I LL t 200 » 10 ) 

EQUIVALENCE UFILLUtl ) tIFILE (4001 ) ) 

EQUIVALENCE (IF ILE < 1 ) *TEMP ( 1 ) ) # ( IF ILE ( 1201 ) *FEMP ( 1 ) ) 

EQUIVAL ENCE ( IFILE< 2401 ) * JTEMP( 1 ) ) * (IF ILE ( 360 1 ) *F1 ( 1 ) ) 

EQUIVALENCE ( IFILE(4B01 )*F0R(1) ) 

DATA F/ * LOAD «*«TIME ***FORCE **'RESP 1 * ’ PLOT 1 * 1 CYCLE */ 
901 FORMAT ( A6 » 3 1 4 ) - 

902 FORMAT ( /// * NO LOAD CARD - QUITO 

903 FORMAT ( /// 1 ERROR - RESPONSES ON TAPES»*2I4) 

904 FORM AT (/ / /* ERROR - MOUNTED T APE 1 * 14* * NEED TAPE * » 14) 

905 FORMAT ( / // * RESPONSES 1 * 14 * 1 T0SI4*' LOADED FROM UNIT’*I3) 

906 FORMAT ( ///* USING FILES WITHOUT LOADING FROM TAPE ‘ ) 

907 FORMAT ( // / 1 ERROR* UNlT" f *I4) 

908 FORMAT ( / // 1 NO TIME CARD - QUITU 

909 FORMAT!///* SET TO RUN'tlSt* POINTS AT’*I3** MS*/ 



I OUT* 6 



Figure A-2 
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WRITE! I OUT • 909 ) N PT $ * N MSt MT I ME 

IF (NPTS.QT.l. AND. NPTS.LT.1201.AND.MTIME.LT. 6001 > GO TO 44 
WRI TE ( I OUT *910 ) 

GO TO _99 

C SET UP FORCING FUNCTIONS 

C READ FORCE CARD 

44 READ ( I IN #911 ) FF » (MT ( I) »DT U > > I ■ I * 6 ) 

IF(FF.EQ.FO) ) GO TO 50 
WRITE ( TOUT *912 ) 

GO TO 99 

30 IF(MT(1UGT*0) GO TO 52 
WRITE ( IOUT *913) 

GO TO 70 : 

52 NSET ®MT ( 1 ) 

MTT=N5ET 


DO 2 I = 1 • N I 


2 

64 

SC(I)*0. 

DO 54 1*1 *MTT 
TFMP! n«I#DT( 1 ) 



DO 56 J»2*6 

I F ( MT ( J ) .LT.l ) GO TO 60 

MTJeMTU) 



mtt«mtt+mtj 

IFIMTT.GT.1200) GO TO 58 
DO 56 I-l.MTJ 




56 

NSET = NSET + 1 

TEMP { NSE T)*TEMP(NSET-1) +DT { J ) 
GO TO j&JJ 




50 WRITE! IOUT *914 ) 
GO TO 99 


C INTERPOLAT E..T1 WES. 

60 DO 604 I»1*NPTS 
IS-I 

_ TG = I »NMS 

1 F (TG.GT « TEMP ( 1 ) ) GO TO 606 
JTEMPU>*1 

604 Fit I) « TG/TEWP.ll 

GO TO 612 

606 J® 2 

: DO 61 0 Lr l-SimtS 

TG»I*NMS 

607 IF ( TG.LE.TEMP < J ) ) GO TO 608 

J g 

GO TO 607 
600 JTEMP ( I >-J 

Fl(n«{TG-TEMP(J-l) ) / 1 TEMP ( Jl-TEMP ( J-l > ) 

610 CONTINUE 
612 CONTINUE 

6 JL_READ..LtiJ4 i$JLSUL t.L*SCL 

IFML.LT.l) GO TO 70 


Figure A-2 ( Cont r d ) 

Program II * Fortran Listings 


sc( i u-«sa — 

I F ( 1 L«LE«N I ) GO TO 62 
WRITER I OUT *916 ) 

— G O TO 93 

62 READ(IIN*917MFEMPtI>*I«l*MTT> 

DO 66 I- 1 * NPT S 

— J«JTEMP( n 

FOR ! I >«FEMP ( J ) #F1 ! I 1 

I F < J« EQ« 1 ) GO TO 66 

FORt I ) -FORI I ) +FEMP1 J-l ) «( U--EJ UJ.1 

66 CONTINUE 

WRITS 13 ML) 1 FORI I ) tl-ltNPTS) 

WRITE! IQUT *929 ) I LiNMSiNMSt! FORI! ) » 1*1 tNP . T S i 

GO TO 61 

70 CONTINUE 

-C FINISHED READING^ INTERPOLATING AND STORING FORCE R 

C REAO RESP CARD 

71 CONTINUE 

READ ( IIN * 918 t FF > IRA » I Rfl i.llf l.L. f . FR Efl. — 

IF ( FF *NE • F! 4 } ) GO TO 99 

IF(IRA.GE.IRX.AND*IRB*LE.lRV#AND*IRAiLE.IRB) GO TO 72 

WR I T E < IQ UT 1 9 1 9 ) I RA > I RB , 

GO TO 99 

72 WRITE! IOUT*9190> IRA * I RB 

1F(I2FIL>GT*0«AND»1 2F It • LL*A _ ) GO TO 74 

WRITE! IQUT * 920 ) 

GO TO 99 

74 GO TO t 76 * 78 » 8 0 ) * 1 2F I L 

76 WRITE!I0UT*921> 

GO TO 82 

78 WR I TE t TOUT >922 ) 

GO TO 82 

80 WRITE! IOUT *923 ) FREQ 

82. CONTINUE : 

C READ PLOT PRINT CARD 

83 READ ! 1 1 N *927 ) FF*NPL 

I F ! FF ■ EQ * F ! 5 ) ) GO TO 84 

WRITE! I OUT *9 30 ) 

GO TO 99 

84 CONTIN UE 

DO 8.00 I R» I RA * I RB 

DO 110 I * 1 * NP T S 

no suMtn -o* 

SUM ! NPTS+1 ) -0 • 

IRZ-IR-((IR-1)/NR)*NR 

READ 1 1 f 1 RZ ) STAT 

DO 700 IL»1*NI 
IF!SC!IL)«EQ*0») GO TO 700 

I Fa U I RZ-l) »Nl+( IL".l.n *LQ : 

DO 120 1-1,10 
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l£E=XF+I — 

120 READ <2* IFF) (I F ILL < J 1 1 > * J«1 ♦ 200 ) 

MPAC=STAT (3* IL)+«1 

_ I F ( I Z F T L_ i_E Q » 3 ) G O TP..3L22_ 

IF(MTIME.LT*MPAC> MPAC-MTIME 
122 J»0 

— DO 124 , Ij »_1*MEAC 

K®4001+ ! 1-1) /3 

IFILE! I > “FLO (J.*12#IFILE<K))-2048 

J_= J+.L2 

IF ( J*EO» 36 ) J a 0 
124 CONTINUE 

0 ) t I Z F I L 

130 ZER*STAT(2tIU*NMS 
GO TO 132 

131 ZFRsSTAT(4»IL.)»NM5 

132 MPTS°NPT S 

I F ! MPAC«LT #MT IME J MPTS»MPAC /NMS 

GO TO 113 i 

150 CONTINUE 
MPTS=NPTS 

IFtMPAC « LT«MTIME. L MPT5*MPAC/NMS 

NEED*MPTS*NMS 

CALL FILTC IFILEfMPACiNEED#FREQ) 

lERmSU 

133 CONTINUE 

DO 140 I=1*MPTS 

I T F ( T UP 

DO 138 J«1*NMS 
K«n-1)#NMS+J 

1 3 fi I TF M _L= I TF tl ) + 1 F I LE.(Xi 

140 YITFU )-ITF(I J-ZER 
400 CONTINUE 

RFAD ( 3 1 I L ) (FOR ( I )_ »_l_ a 1 tMPT-SJ 

MPTSS«MPTS+1 
DO 402 WtMPTSS 

402 RFSP ( I )«Q> 

DO 404 I = 1 * MPT S 
JM=MPTSS-I 

DO 404 L*1j^ — 

jp * J + 1 

404 RESP(JP) a RESP(JP)+FOR( I ) *Y I TF { J ) 

S.C AL£.gATAl i-LiJLLyuA01»SC( IL1 

DO 406 I »1 iMPTSS 
RESP ( I )=RE$P( I )*SCALE 

406 SUM( I )« SU M|J_U -_R£j P( I ) 

WRITE(IOUT#925) IR* IL*NMS* <RESP( II * I*2*MPTSS) 
IF {NPL*EQ*2 ) CALL PLOTB ( IRt IL#MPTSS*NMS*RESP) 

7 00 CONTINUE 

WRITE! I OUT *9 26 ) ( SUM l I ) » I*2*MPTSS> 
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.♦mstSUMJ 


800 CONTINUE 
READ RECYCLE CARD 

DP AH IT T N I PP*NCY 


IF<FF.NE«FI6) ) GO TO 99 
WRITE (IOUT, 9 2 8) 

r,0 TO ( 1 »3&iA4.i_61j,7 1 iSS. Li-NCY . 

99 CONTINUE 
CALL PLOT C 



SUBROUTINE PLOTS ( I R * IL #MPT »NMS * Y ) 
DIMENSION Y ( 1 ) 






RETURN 
ENTRY PLOTA 

PLOT IN IT IALIZATION HERE 
RETURN 
entry PLOTC 

PLOT TERM INATION HERE 

RETURN 

END 
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kEOHUS fcAS-TX»EASU- 



SUBROUTI NE FASTI < X *Y *N * Z *W *S ) TW002260 

FAST FOURIER TRANSFORM OF COMPLEX DATA BY DAVID IVES* GRUMMAN TW002270 

INPUT VALUES ( COMPLEX ) iwoozzsil 

Y**«.*N OUTPUT VALUES (COMPLEX) TWO0229Q 

#•••• .NUMBER OF VALUES TWOG230Q 

Z ui • aDUMMY STORAGE OF LENGTH 2N (COMPLEX) IkOfi2_31iL 

W* • * # # DUMMY STORAGE OF LENGTH N (COMPLEX) TW002320 

5« • • * #SI GN CONTROLLING DIRECTION OF TRANSFORM TW002330 

...IHIS PRODUCES * OUTPUT Y» FROM » INPUT X»t WHERE TW002 340 

«*####*#####*###**###***#*###*###############*# ###*##########*####7^002 350 
K*N TW002360 

H J ) *5UM X tK) »EXP ( S I GNi 1 i » S ) » I » 2«P 1 » ( J-l ) » ( K- 1 ) / N ) TW00 2370„ 

K*1 TW002380 

* ###****# ####**#* *### ** *-*■#**** ****** ***** ** * ** ******************** f ^ o 0 2 3 9 0 

J^ITH I a SQRT ( *1) * S» + l» OR S*-l«» AND Pl=3»14159 TWQQ2400 

COMPLEX NUMBERS ARE HANDLED IN FORTRAN 4 CONVENTION* NAMELY THE TW002410 
REAL AND IMAGINARY PARTS ARE STORED IN ALTERNATE CELLS* STARTING TW002420 

MITH T HE R EAL . P A RT OF .. T HE. F IRST LQC A T LQIitJEIila MQOMlfiL. 

DIMENSION X ( 1 ) * 2 ( 1 ) *W ( 1 ) * Y ( 1 ) TW002440 


MOD ( J*K ) “J-( J/K ) #K TW0Q2450 

DO 1 I a 1 *N TW0Q246Q 

W ( 2*1-1 ) «C0$ ( ( 6.2831853071 7959/N )*< 1-1 1 ) TW002470 

W l 2* 1 1 * SIGN(1. *S)*SIN( (6*28318530717959/N)*< 1-1) ) TW002480 


Z I2»I-l>«Xt2«I- l) 

1 2 ( 2# I ) *X < 2#I ) 

I D“N 

DO 4 J»1*N 

IF(ID-l) 5*5*15 
15 CONTINUE 

DO 2 IX«2*I P ; 

I F ( M0D( I D » I X ) ) 3*3*2 

2 CONTINUE 

3 XP*1.P /I X 

IS-N/ID 

DO 4 Ll*l * I S 


TW002490 

TW002500 

TW002510 

TWOQ252Q 

TW002530 

TW002540 

TW002550 

TW002560 

TW002570 

TW002580 

TW002590 

TW002600 


D O 4 L *1 1 I D TW00261Q 



JM«<M0D(L+{L1-1)*ID*IX*N)+M0D< J+l*2 >*N>*2 
JP«(L+<Ll-l )*ID+M0D< Jt2)*N)*2 
Z (JP-1 )*Z ( JM-1 ) 

TW002620 

TWOQ263Q 

TW002640 


2UP)-Z( JM) 

TW002650 


DO 4 IH*2 * I X 

TW002660 


IG» ( MOD ( (L1-1)#ID*( IH-1 ) *N ) +1 ) *2 

TW002670 


lU«JM+{IH-l)*iD#2 

TW002680 


2 ( JP-1 > ■ Z < JP- 1 ) +Z ( I U-l ) #W ( I G-l > -Z ( I U > *W ( I G > 

TW002690 

4. 


TW002700 

5 

DO 6 1*1 *N 

TW002710 


K* 2#(MOD( J+1*2)#N+I )-l 



Y(2*I-1)-Z(K) : 

TW002730 


k*k+i 

TW002740 

6 

Y(2*n*zno 

TW002750 


RETURN 

TW002760 


END 
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♦ (RUN CARO) 
*ASG .T 1 * F2/S7/4 

♦ ASG * T 2 * F2 /7/P OS/7 
*ASG*^TB 1 0 *8C . A065 1 5 

♦ A S G T B ~ U ^SC rrARES 

♦MOVE 11., 1 

♦ASG.TB 12, 8C * T APE6 

♦ MOVE 12 * * 1 

♦ ASG * TB 20 * 8C * SCR ATCH 


FORTRAN DECKS 


♦MAP, I .MAP, .PROG 
SEG TOP* 

IN ITF1 "T 

LIB MSC*LOCALIB 
♦ XQT .PROG 

lo ad "nsr i “ ■ ro “ 

STATUS i 10 

RUN 11 5 1 10 S 

j2- 500V 1X43 ; 63 g # - I243 - 624 ." TO 92V I TTO. 653. 64 (TV* 1 79. 1 79 . " 

RUN 12 6 1 10 5 

12500. 1343. 750. 1460. 733. 1284. 1305. 767. 752. 210. 210. 

PLOT’ 1 10 5 6 500 4 2 O " 

STATUS 1 10 

OUMP 20 1 10 


*PMD. E 
♦ FIN 


* = 7/8 PUNCH 


Figure B-l 
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RECORD NO* *69 


LO AO NO. 


IMPULSE TRACK 


14 4 -S 46 18- 


RESPONSE NO 


16 17 19 21 22 

Id -19 22- -22 -2-1 

16 14 14 1 5 IS 

iq — -2*. — is — 42 


a lo 

|2 H 

S 3 

- 9 a- 


21 17 15 13 


-80— -K 2 A 78 857 6*4 - *92 — 487 h82 

9 


9 9 

40 1 0 


9 B 9 

44 LS — — 1_4. 

12 12 12 

14] - -a a - 

18 17 II 

4-y 14 14 

9 12 19 


1 1 
I 2 

17 
1 5 

1 1 

| 4 - - 

9 

1 

12 

44_ 

13 

I 2 

1 2 

1 3 

1 3 

44- 

10 

-44- 

13 

-4 4 

IS 

12 

1 I 

9 

1 1 


44 

- 444 

42 IS 


-1 9 11 

— a 4. 

4 2 4 

— S 1- -| 

11 14 13 

-3 2r 1-0 

4 -6 >8 

22 1 3 •! 

5 17 18 

5 7 9 

0 1 2 

-4 - 8 -9 


13 6 

4-0 1 4 

-8 2 

• 1 —4 1 

18 13 

9 A 

2 t 


10 11 
45 42- 

12 13 

13 14- 

16 17 


13 14 

40— — 40- 
1 1 9 

2- 

13 11 

12 1 l - 


3 

.■ 6 - 

6 

14 

L 2 

13 
1 4 

1 1 

u 

10 

10 

s. 

7 

4 

4 

4 

8 

13 

IS 

14 

13 

1 3 

-4- 

1 

I 4 

1 7 

14 

1 l 

111 

a 

0 

□ 

1 1 

12 

12 

I t 

10 

2 

- 2 - 

— 0 — 

-40 - - 

-0 

• — -o_ 

7 

s 

s 
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RESPONSE NO. 2 


RESPONSE NO. 3 


22 I 

-33 4 

2** 1 


t 0 

— 8 

I 1 

-40 

I 3 


1 2 — 4-3 
B 7 


-444 +0~ 

12 12 

IS 1-4 — 


•33 -29 


RESPONSE NO. 


1-2 — 

10 

— I-Z- 
15 

US — 

12 

12 

1 7 

14 - - 
1 A 

—44 14 

-1 -1 

- - 1-0- 

2 

g 

8 

— 1 - 

1 3 

- 4 - 

7 

IS — 

- 75 - 

17 - 

— 1-2 

- -U 

- 1 1 

15 12 

6 

-1 - 

-8 

2 

20 

12 

12 

t u 

8 

9 

a 6 

9 

10 

1 0 

5 

21 _ 

1 ~ 



— 1 1 

42 

12 

- rt- - 8 

1 

0 

2 

- - 8 


- - 30 - 

— — 5 

-5 

- 4 

- 8 

9 

1 1 

10 

1 1 

8 

1 

2 

A 

6 

3 

8 

— 34 - 

_.4 

- 3 — 

1 

— • 3 

..-Si .. 

— 1 - 

— 5 - 

«* 

15 

i a 

15 

l 1 

10 

A 

1 

-38 

10 — 

- 40 - 

-8- 

— 8 — 

11 

- - 9 - - 

-4 2 

-52 

1 * 

1 7 

17 

IS 

13 

30 

31 


RESPONSE NO. 


RESPONSE NO. 8 

4-*3 2 0 4 218 >4? 25-8 

27* 277 272 268 27a 

Z8H 2*4- 390 — -2*2—772- 

2*9 23 1 222 213 200 

4-8-3 3-4-3 34S 127 10* - 

88 64 *7 31 20 


13 13 IS 

20 21 1 * 


-2 -3 


7 12 1* 

3*7 
- 8 8 7- 

13 18 22 

—4 4 *- 

23 20 13 

— 5 * 44 

21 17 IS 

-4-S 1-4- -4-3 

18 18 18 

— S 6 1 - 

U 18 22 


1-2 2- 

I 10 


2 * 33 3 S 

—1 7 4 -* 4-8-- - 

a s 7 


13 

— . 12 
7 

- 7 _ 

1* 


IS -61 -92 -128 

- —48 .2-17 — -225- -740 

-1 -263 -285 -304 

7 -333 — 2**-~3*7 

22 -31* -295 -273 

9 -. 23 8. —24 -4 -4 7 2 - 

13 - 1 0 e> -103 -100 

8 -* 5 23 

4 144 1-44 1 7 2- 

12 219 235 252 


* 300 29a 270 

2 3 22* 224 — 71-7 

IS 137 118 105 

5* 40 •74-- 

II -* -17 -35 

t 1 -28[1 -237 -237 

|_Z— _-2 9*_-22G --32 3 

7 -300 -311 -325 

- -13 -2 X3 - 24 5— • 2 28 

13 -20S -175 -139 

24 -42- — —9-8 — — 4-00— 

37 -41 -83 -26 

46 1-8 — 48- 17 -I 

6 7 7 7 

— 4 7 1 4- 4 1- 


-|44 -198 
-274- -244 
-315 -313 
-34-8 -322 - 
-257 -289 

n i 4 i I j \ i „ m 

— » - ■ X % r 

-90 -77 

—22 —21 - 
87 76 

1*0 20* 
263 270 

- 2 8 7 7*4 

2S0 235 

-- 44 S- -- 4 - 6 S- - 
9* 78 

_ _i 3 5_ 

-SA -83 

-314 —22-8- 
-2*9 -271 
-“24 7 — -20-1 - 
-32S —305 
-223 — 249 — 
-107 -92 

--—9.4 — -79 
-10 6 

8 - IS 14 

6 709 -3 

2 39 - 22 


RESPONSE NO. 
-4 7 14 |S — 

25 29 29 

-1-5- IS — 29— 

26 22 17 

-29- - 28 25 - 

26 26 25 

-2 4 28 22— 

21 20 20 

22 31 32 - 

18 19 28 


78- 27 

29 27 


20 20 
31 32 

1 9 28 

22 70 

27 21 

-30- — -34- 

20 20 


1* 1 a 20- - 

25 28 28 

27 26 5-3 — 

16 17 21 

. — 2 -S 23 3* — 

20 25 27 

- 4-8- 14 1-3- - 


24 34- — 24- 


l 7 — | *1 


28 

27 

26 

15 

-18 

23- 

28 

18 

15 

gy 

-25- 

1 4 

-29 

-15 

19 

-43 - 

19 

S- 

23 

29 

37 


Figure B-2 (Co 
Program I - Record Out 


£ 





( 



IMPULSE AREA IN LB-SEO 1.6RR965 
4)UR AX44)N-^. -fU-SPaALS£^ 4-N -SECON-OS* -430 2* - 


S3 mi A 4 3X- O F R EC OR D A V E R AS5S “ 

___&££-• -4-M8 _R4 R-2 63 6ft R5 R6 - *7 R® R9 *1° 


120 ID V |ft 7 2 10 2 10 15 22 27 

1. 30... 2 40 43 7 1 1-3 5 $ 16 — 23 ~X5- 

1 RO 13 11 13 5 1 9 ft 1 1 1 5 23 26 

1 5 q 7 j_j. 4 — — 1 14 - 2 - -- 9 1 R 22 2 7 

160 1 S86 1616 1583 161R 1617 1596 176R 1601 1615 1793 1802 

373 4586 1.414 -1 -583 184ft — 464 7 4 546 17 6-R 460 1 - 16 45 47 94 4 -802 

180 1585 1617 1583 161R 1617 1595 1766 1601 1615 1793 1802 

W 45-8-6 1-64-7 1 -5 - 8 3 4-64-4 1 6 17 — *5-9-5 4767 4-6454 — -46-4 5 479 3 4802- 

200 1587 1617 1583 161H 1617 1596 1768 1601 1615 1793 1802 

2 1 q—4 566 4647 4563 4644 — *64 7 1695 - 4767 460 1 464 5 1793 I 802 

220 1587 1617 1583 161R 1617 1595 1767 1601 1615 179R 1802 

240 456-6 4647 4563 46 44- — 4647— 4 59-4 4X68 4601 1 61 5 I 7 93 4503 

2 RO 1 586 16 17 1 583 161R 161 7 1595 1768 1601 1615 1793 1802 

350 4567 4647 4563 4644 4647 4595 4X66 4662 464-5— 4763 4803- 

260 1586 1617 1583 161R 1617 1595 1766 1602 1615 1793 IBQ2 

47Q. — 4&-fr4 464X — 16 8 3 —4644 164 7 4594- 4767- -4604 4 64 5 4 794 180 2 

280 1587 1617 1583 1614 1617 1595 1769 1601 1616 179R 1802 

. -2-90 1567 4647 1 6 8 3 46-14- - 46 47 15 9 ft 4? 69- 160 2 1615 i 7 9 R 18 0 2 


300 1 586 161 7 1583 1 61 R 1617 1595 1 7 7 1 1 602 161 6 179R 1803 

34_0 43 8 -9 ft 6- 49 47 44 - - 9 - * 5 — 12 

320 13 8 9 8 6 1ft 16 10 9 1 8 1 1 

3 30 9 ft -9 6 6 -l-ft — 46— 40 — 9 1 6 14 

3R0 | R 8 9 8 6 1ft 18 10 9 |8 11 

35-0- — 8 8 40 6 6 — 43 16 4 0 9 18 44 

*♦61 1ft 7 9 8 6 1R 20 10 9 1 7 40 

R32 42 7 9 7- 6- 43 2-0- 4-0 X - 48 -4-0- 


R63 R87 513 R90 503 R61 R99 R91 R67 R39 R95 519 

ft 4 ft m 6 --48 6 0 2 3- - 40- -15 2-1 - 22 

R65 10 9 16 1 10 ft 11 -19 21 28 

4 44 8- 4-1 44 7- 2- 4 0— — 7- ■■ - 9 - -2ft 23 — 2-7- 

ft 47 1ft 9 16 6 Q 11 7 9 -26 22 25 

j^4B 9 1_0 | 6 7 — j. -37 6 — 1-0- — --—-23 2-2 — 25- 

R69 15 9 9 6 3 1 1 ft 12 -19 22 25 

569 3 -X 46 7 -4 42 6 9 4 9 - -22 25 

570 0 12 7 5 8 13 5 7 15 22 2ft 

5X4 Q 40 -46 7 43— 41 3 X 46 21 26 

572 0 9 18 7-21 2 5 7 16 21 27 

5 7 3 0 0 0 0 0 0 Q - -Q- O 0 0- 


FUNCTJON CARD - STATUS 1 10 0 00 0 0 00 0 


Figure B-3 

Program I - Record Stannary 




Pigurs R-4 

Program I - ITF Listing 


s 

* 

s 





# (RUN CARD) 

#ASG»T i*F2/2//4 


TvtcTl mrm . 1 JrTmWk u*UWMi 


#A$G»T 3#F2/30//30 
#ASG*TB 11 * 0C *A065 15 



# MAP # I #MAP**PROG 
* 


IN RESP 

LIB MSC*LOCALIB 
»XQT •PROG 


DATA DECK. 


*PHD 

#FIN 



Figure I 
Program IX - 1 






2 >U 




1 l& 

2 . , • 

63 2. > 'I' 1000 






6 20 25 25 59 234 695 1056 1189 HOC 920 698 813 1113 1430 1610 

1789 2044 2375 2690 2963 3109 3178 3261 3355 3524 3756 3925 3921 3886 3899 3914 

3898 3887 3896 3891 3876 3896 3908 3892 3909 3918 3928 3895 3880 3884 3894 3904 

3920 3920 3920 3920 3920 3920 3920 3920 3920 3920 3920 3920 3920 3920 3920 

1* 

00 0 0 34 235 567 992 1037 1013 1025 1048 1094 1206 1341 1517 

1750 2091 2354 2435 2513 2603 2701 2814 2926 2987 3067 3183 3243 3264 3274 3289 

3201 3248 3222 3207 3248 3259 3235 3168 3185 3202 3210 3224 3225 3221 3203 3231 

3233 3235 3235 3235 3239 3235 3235 3235 3235 3235 3235 3235 3235 3235 3235 3235 

6 . «!• 

0 0 0 0 -7 -70 -243 -876-1271-1312-1224-1262-1339-1483-1678-1795 

-1920-2081-2293-2455-2563-2624-2703-2831-2976—3167-3350-3464-3361-3315-3360-3423 
-3442-3403-3336-3217-3267-3292-3294-3290-3270-3265-3268-3269-3251-3234-3230-3238 
—3243— 3244— 3244— 3244— 3244— 3244-3244— 3244-3244— 3244—3244— 32^4-3244-3244— 3244— 3244 


5 !• 


0 Vo 0 0 0 -32 -224 -746-1085-1319-1369 t1209-U64-1280~1524-1824 
—2049— 2210— 2302— 2497— 2773— 3038-3193-3210-3324— 3504— 36T4— 3709— 3602*^3513— 3586— 3614 
-3626-3633-3634-3634-3611-3628-3645-3631-3629-3625-3608-3595-3624-3637-3606-3581 

- 3580-35 87-3589-3589-3589-35B9-3589-35B9-35 *9-35*9-3589-35 89-35 39-3 589-35 9 9-35 89 


8 1 . 

0000 -94 -298 -779-1161-1299-1332-1263-1155-1270-1537-1813-2069 

- 221 4-23H-2494-2766-2999-3167-328J-33 17-3432-3548-3621 -3? 72-3*8 3^345^^344 8-3462 
: -#40-f^4**^46 , 4-3464-3464-3464-3464-3 464- 3464-3464-3464-3 464-3464-3464-346 4-3464 
<*3464— 3464— 3464*^3464— 3464— 3464— 3464— 3464-3464— 3464— 3464— 3464-3464— 3464-3464-3464 
9 1* 

0 0 0 0 48 218 484 838 1060 1130 1176 1230 1176 1250 1386 1571 

1823 2123 2432 2717 2828 2885 2940 3137 3282 3422 334C 3296 3438 3416 3370 3361 

-2342 3394 3479 3424 3393 3368 3348 3381 3385 3385 3385 3385 3365 3365 3385 3385 

3*89 3382 338 5 338? 3385 3365 3389 3385 3385 3385 3385 3385 3383 3385 3385 3385 

•ml:. ; : : , - 

0 0 0 0 30 283 834 1272 1039 943 1170 1422 1415 1418 1538 1659 

1879 2319 2605 2605 2578 2679 2866 3039 3211 3353 3473 3629 3564 3471 3534 3560 

3360...M.36.. 3 543 3 5 7 6 . 3575 3 ,569. 3 556 ,. 3 .5 65 3575 3555 353 . 7 3 5 49,. 3 5 3 , 8 . 35 6 Q 3.5,5j>„ 

3546 3545 3545 3545 3545 3545 3545 3545 3545 3545 3545 3545 3545 3545 3545 3545 

ii u 

1671 4937 T213 9045 7664 5202 4472 4691 6050 7506 7758 7210 7295 7308 8358 8659 

8690 8880 8992 9107 9750 9967 9918100771076611350115051118811047109701090910757 

10781109001081610579105301064710873110751123811252111101101610956108871077310682 
106631 0730 107591075 91Q759 1Q759107 59 10759 10759 107 59107591 07591 075 9 1Q7 59 1075 9 10759 

v" "f - 7 -• 

-166 lltl 4277 6339 7737 773 2 6 044 515 1 5374 6585 8255 8897 8 1 35 796 1 85 0 4 9 0 94 
9324 9624 986510108103601056910837112941152811457113041128711413115661177811651 
114331113111061112111138311465114381152611662X1772118161165911378111781111711283 
11322113221131911319113191131911319113191131911319113191131911319113191131911319 


_JL3 l»_ 

367 2335 6370 7997 7999 6662 4957 4990 6883 8338 8068 6831 7026 8305 9093 9854 
8457 8165 8693 99231037210286100501012910836115511 17971 15031;l240ll587il9341 1754 
11247108731102411348115431156211474114281152211593115391133311227112311130311276 
11 1571 117811200 1 120011200 11200 112001 1200 1 1200 11200 11200 11200 11200 11 200 112001 1200 

14 1. 

660 3905 7051 8381 8121 6555 5149 4792 6076 7502 6640 82B9 7467 7727 9024 9393 
T3393 943$ 9507 97211020010414106341 1239 11429 U358 l 1163 1 0987111401 1318 1131 81 1274 

, v -V- ; < • ' - 7 ; • ‘ 



( 

( 

( 

< 

< 

( 

< 

( 

( 

( 

< 

( 

( 

( 

"3 

X 

l 

k 



;■ Figure B-6 

j ” • *• v 

Program II - Beta Deck 


r 

1 



— 


>** 



114951149511495114951149511495114961149511495 

■ " ' " L - ‘ 1 “ T M * ^ " ™- nr 7 i'! *■’ - * * • » '< 


. • ■ - : V. - ■ - ■ -r • .y . r • . • ... O;-;. "• -a . • 

*76 014710691 9357 7709.6220 6610 *61710650 9076 7 50 2 807**3*14-031 ?l0 550 
9428 9253 9627114 1 41 2006 1 1406 10081 1 0 1 19 1 175? 125 69 123 1 51 1 40 411378 1 18071208 01 1 9Z3 


my 
^‘VN '.-‘i 

1 •' ( 


11611113601136811808118621142410913109931141311413111921103411029111341119111115 
1 104810995 10980109801098010900109801 0980 10980 1099 0109801 0980 10980109801098 01 0900 
20 1. 


257 1275 3901 906710490 9970 7768 5810 7431 98 12 11 12010336 9385 8866 8980 9717 
1G751U2341 1271111521 1163 113381 17&012238 12544123941 18651132711495119691219412 188 
1 2Q1611772 116351187312047 119501:164911576 11727 118621 19361 1874117481 172 11 176411775 


11 72 3116221158511585115851158511585115851158511585115851158511585115851158511585 
21 1 

1088 4353 909310810 9930 7881 5924 6394 B4QP 931410325 9645 6322 7729 7984 9295 


,*68010088102631082511 13610070103601025 11114e 118621 19331 1380 11036108951X05 011527 
li¥97ii778 115401 1245 109921091 0110591 1296113241 1199110541 1012 1109011201 1131 71 1291 
11016108771085010050108501085010850108501085010850108501085010850108501085010850 


22 !• 

558 2791 6812 997810348 8733 6245 5379 74921016910738 9567 8690 8402 8552 9174 
101491039010106 9999 10333 107171 1048 1 U2 1 1 16261 1477j IQ IQ 1 0567 105 83 1 1 002 1 155 61 17 18 


rrl|4noai910766 110781130011261 11008109141093211014 11 1591128311065109231093211060 
1345 11345 11 345 11345 11 345 11 345 11345 11345 11 345 11 345 11 345113451 1*4 511345 

_24 Li ' ■ . • -l, 


0000 -26 -183 -360 -634-1003-1239-1288-1253-1238-1279-1311-1363 
-1548-1789-2045-2268-2442-2542-2601-2676-2864-3060-3238-3416-3465-3426-3444-3503 
-3555-3590-3608-3578-3584-3598-3611-3605-3575-3556-3542-3544-3552-3534-3517-3514 


*3320-3340-3544 -3544-35 44-3544-33 44-3 5,44-3 54 4-3544 -3544-3544-3 ^^^^^3^^4^8544 
-20- 24 5 .69 2 1 QS 8 . 1298 1 3Q7 12 24 11 $4 l gp9 13 34 , 1 945 . 1 7 77 , 20 1 2 


2207 2409 2708 2904 3087 3217 3272 3298 3331 3364 3406 3451 3415 3371 3410 3465 
3456 3492 3561 3540 3454 3419 3422 3425 3414 3398 3418 3474 3537 3582 3584 3505 
3431 3404 3400 3400 3400 3400 3400 3400 3400 3400 3400 3400 3400 3400 3400 3400 


06- 


'“ 1 * • -■ ■ • •’ ‘' V '-V ■■ 

* < *4 -12 -122 -463 -890-1219—1398-1427-1343-1215-12^414*6-1716*1911—2051 

2 ^ 1r248l-2Sa7-L636*27l 4 *Lfl3:7-l Ml* 3.L ^ j^289 - J. 34 .6-33J3 1-3262-3^ f!?3g G:i .r. a?l .lra2 7 .0 


-3229-3202-3196-3217-3227-3234-3241-3275-3235-3208-3100-3174-3159-3150-3117-3105 

-3127-3095-3079-3079-3079-3079-3079-3079-3079-3079-3079-3079-3079-3079-3079-3079 

29 


0 000-7 -79-297-865-1259-1376-1298-1265-1349-1536-1853-2093 
-1006-2407-2642-2812-2906-2979-3 164-3396-3483-3539-35 29-35:62-3352-0261*3328-3350 
^ly^ 34 3 3r3 4 S9 r 34 .51 - 34Qa- 33 89- 34 O3-3 4O O * 33?5-3394-3392-33,72-3098*S ^98 , *33, 94-3412 


-3411-3378-3364-3364-3364-3364-3364-3364-3364-3364-3364-3364-3364-3364-3364-3364 
35 1« 

-3 -4 -4 -4 59 337 968 1262 1310 1256 1220 1158 1304 1508 1734 1945 


.0147V 02.95 2504 2697 2758 2824 2911 2970 3147 3280 3250 3105 3^^143 0198 3032 
*0*7 3145 3050 3129 3187 3195 3243 3278 3189 3153 3147 3159 3131 0098 3091 3064 
3047 3055 3060 3060 3060 3060 3060 3060 3060 3060 3060 3060 3060 3060 3060 3060 

36 -1* 

0 0 0 0 6 65 269 912 1277 1418 1408 1298 1213 1328 1563 1815 
1967 2161 2387 2662 2908 3074 3231 3336 3432 3655 3771 3700 3668 3659 3620 360fr 
J3033 3689 3752 3805 3832 3821 3766 3720773799 3774 3718 3812 3809 3788 3801 3868 
3787 3705 3795 3795 3795 3795 3795 3795 3795 3795 3795 3795 379$, 3795 *799. 3795 


H 


!( 







































FIGURE C-2 (CONT'D) 
PROGRAM I - SUBROUTINE STORED 

















FIGURE C-4 

PROGRAM II - MAIN PROGRAM 
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FIGURE C-4 (CONT'D) 
PROGRAM II - MAIN PROGRAM 











FIGURE C-5 

PROGRAM II - SUBROUTINE FIUT 





























